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(C)  Aii  exiierinientnl  study  uf  coherent,  matched  ifitortoelTrrtrpios  for suhmarine-to-siibinnrino 
communication  was  conducted.  Coherent,  matched  filter  techniques  offer  three  advantages 
over  existing  systems.  First,  coherent  integration  allows  useful  system  operation  at  signal-to- 
noise  ratios  significantly  below  0  dl>.  Incoherent  systems  experience  a  threshold  at  a  0-dl> 
signal-to-noise  ratio,  which' causes  system  performance  to  deteriorate  rapidly.  Second,  filter¬ 
ing  of  tin*  received  information  symbol  with  a  filter  matched  to  the*  received  symbol  wave¬ 
form  reduces  the-  effects  of  both  noise  and  multipath.  Finally,  those  techniques  are  compatible 
with  a  wideband,  randomized  transmission  format,  which  reduces  tin*  detectability  of  the 
signal  by  unintended  receivers.  The  research  program  was  conducted  in  two  stages.  The  first 
stage  evaluated  communication  system  performance  over  7-  and  42-iuui  fixed-site  ranges. 
Over  4,000  hours  of  experimental  data  were  obtained  and  analyzed  to  ensure  the  statistical 
significance  of  the  measurements.  Ke<u1t~  from  those  data  showed  that  reliable  communica¬ 
tion  could  lx  obtained  with  signal-to-noise  ratios  below  0  db.  A  typical  cx|K*rinicntal  system 
(Mf»B)  transmitted  0.f»2.">  bit/sec  in  a  100-hz  band  centered  •"»  420  hz.  Over  a  42-nmi  patl\ 
this  system  yielded  a  bit  error  probability  of  0.01  at  an  input  signal-to-noise  ratio  of  —0  db. 
Such  pert  o  nuance  is  within  0  db  of  that  obtainable  with  the  optimum  receiver  operating 
through  linear  time  invariant  channel  with  added  white,  (iaussian  noise.  To  determine 
the  applicability  of  the  fixed-site  results  to  the  submarine-to-submarine  communication 
problem,  tin-  second  stage  of  the  program  investigated  the  space-time  stability  of  the  acoustic 
medium  with  a  towed  source.  An  experiment  was  conducted  in  the  deep  ocean  off  Klcnthera, 
B.I.  at  ranges  from  0  to  400  umi.  This  experiment  indicated  only  a  10  percent  decorrelation 
in  the  channel  from  one  00-sec  interval  to  the  next.  Consequently,  tin*  coherent,  matched 
filter  techniques  evaluated  over  the  fixed-site  ranges  are  applicable  to  practical  problems 
involving  moving  platforms.  Based  on  the  results  of  the  experimental  program,  development 
of  a  practical  submarine-to-submarine  communication  system  using  coherent  matched  filter 
techniques  is  proposed.  A  particular  implementation  (M7)  incorporating  a  randomized, 
burst-type  transmission  for  detection  resistance  is  suggested  as  a  basis  for  future  work. 


I.  INTRODUCTION 

(C)  Underwater  acoustic  communication  systems 
based  on  coherent,  matched-filter  techniques  have 
been  studied  experimentally.  Although  the  specific 
goal  of  the  research  was  tactical  subinarine-to- 
subinarine  communication  systems,  the  results  are 
applicable  to  other  situations.  This  paper  describes 
extensive  fixed-site  communication  experiments 
and  an  important  transmission  measurement  that 
encourage  the  immediate  consideration  of  these 
techniques  for  submarine  communication  systems. 

(C)  Coherent,  matched-filter  techniques  offer  three 
distiuct  advantages  for  submarine  communication 
systems.  First,  coherent  integration  allows  satis¬ 
factory  operation  at  low  signal-to-noise  ratios, 


without  the  threshold  effect  that  «  common  to 
incoherent  systems.  Second,  matched-filter  opera¬ 
tion  reduces  the  effects  of  both  noise  and  multipath 
under  varying  propagation  conditions.  And  finally, 
these  techniques  are  compatible  with  a  random¬ 
ized  transmission  format  that  reduces  the  detect¬ 
ability  of  the  communication  signal. 

(C)  The  processing  techniques  employed  in  the 
research  are  based  on  well-established  theory.  The 
acoustic  medium  is  assumed  to  approximate  a 
linear,  time-invariant  channel  with  added  white 
Gaussian  noise.  For  such  a  channel,  the  optimum 
receiver  is  composed  of  a  filter  matched  to  the 
received  signal,  followed  by  n  threshold  device.  The 
matched  filter  can  also  be  shown  to  reduce  inter- 
symbol  interference  caused  by  multipath.  One  of 


CONFIDENTIAL 


I 


CONFIDENTIAL 
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am  sequence  coefficient. 

«(*.  i) 

symbol  composition  coefficient  for  the.*>ith 
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probe  component  of  6th  frame. 

digit  of  the  jth  symbol. 

0|(*.  0 

information  component  of  6th  frame. 

r' 

offset  between  signal  and  rrr  spectral  line*. 

Op(») 

integrated  probe  component. 

r 

integer  offset  between  signal  und  rrr  spectral 

6(0 

biphase  sequence  coefficient. 

lilies. 

r 

average  speed  of  sound. 

r* 

value  of  r  which  maximises  8NR„0I. 

A 

r 

estimated  average  speed  of  sound. 

cross  correlation  between  ar(/)  and  y (1). 

d(k,  0 

decision  on  the  ith  symbol  of  the  6th  frame. 

^  mult 

multiplication  rate. 

d'(k,  0 

preliminary  decision  on  the  ith  symbol  of  the 

«,<0 

signal  component  wuveform. 

4th  frame. 

»(»),  «(4,  i) 

signal  component  waveform  sumple  in  6th 

D 

output  signal-to-noise  ratio. 

interval. 

Hi) 

ith  information  bit. 

5(6,  i) 

discrete  Fourier  transform  of  *(6,  i). 

fc'.- 

error  count  on  probe  component. 

SP,  8P(6) 

signal  power  in  the  6th  measurement  interval. 

lit 

error  count  on  information  component. 

HNR,. 

input  signal-to-noise  ratio. 

K 

cumulative  error  count. 

SNRou, 

output  signal-to-noise  ratio  on  op(tj  vector. 

K{x) 

expected  value  of  I. 

T 

digit  duration. 

/. 

carrier  frequency. 

T„ 

duration  of  probe  frame  plus  information 

h 

Doppler-shifted  carrier  frequency. 

frame. 

ftlk 

receiver  clock  frequency. 

Tlm 

frame  duration. 

ft  <t 

receiver  center  frequency. 

Tm, 

duration  of  measurement  interval. 

fU) 

transmitted  symlsit  waveform. 

T» 

duration  of  segment  of  input. 

til.-,  1) 

transmitted  waveform  for  the  6th  symbol. 

T'„ 

Doppler-distorted  segment  duration, 

6(0 

received  symbol  waveform. 

symbol  duration. 

6(4,  0 

received  waveform  for  the  4th  symbol. 

sequence-period  duration. 

Af  p 

AID  sample  compression  factor. 

u(l) 

transmitted  digit  waveform, 

A'i. 

nurulier  of  sequence  periods  per  frame. 

t/(  0 

unit  step  function. 

a;,. 

numlicr  of  information  symbols  per  frame. 

V 

source  Doppler  velocity. 

Apii 

nuinlier  of  probe  frames  per  measurement 

t»* 

maximum  expected  source  Doppler  velocity. 

interval. 

W 

system  bandwidth. 

A'. 

nuinlier  of  digits  per  symbol. 

W,  p 

low  puss  system  bandwidth. 

A'-p 

numlicr  of  digits  per  sequence  period. 

*,U) 

received  wuveform. 

»'(*,;) 

matched  filter  output  'or  the  jth  digit  of  the 

*'(i) 

sample  of  received  waveform. 

6th  frame. 

x(i),  r(k,i) 

compressed  sumple  of  received  waveform. 

«*, » 

matched  filter  output  for  the  jth  symbol  of 

zAt) 

probe  component  of  received  waveform. 

the  4th  frame. 

*p”(0 

received  waveform  tri  m  infinitely  long  probe. 

w(0 

transmitted  waveform. 

*pE(0 

periodically  extended  probe  component  of 

m,. U) 

prolie  component  of  transmitted  waveform. 

received  waveform. 

vi, ik,  t) 

information  component  of  transmitted  wave¬ 

*i(0 

information  component  of  received  waveform. 

form. 

*(6,  tj 

complex,  compressed  sample  of  received 

v,  U) 

noise  component  waveform. 

waveform. 

mo 

noise  waveform  sample. 

Z(4,  i) 

discrete  Fourier  transform  of  z(6,  i). 

«(o 

noise  average  from  probe  component. 

noise  variance. 

«(i,  6) 

noise  component  of  matched  filter  output. 

,  A* 

chi-squared  probability  distribution- 

A  dp 

degree  of  linear  maximal  sequence. 

'*(«) 

normal  cumulative  probability  distribution. 

JV„ 

number  of  symbols  per  frame. 

y  (4,  i) 

decision  threshold  on  ith  symbol  of  6th  frame. 

Nlm 

number  of  digits  per  frame. 

performance  differential. 

H»tih 

number  of  multiplications  per  frame. 

krtrix) 

frequency  offset  between  signal  and  fft  lines. 

(V, 

number  of  samples  per  sequence  period. 

g(r) 

normalized  cross  correlation  function. 

NP,  NP(6) 

noise  power  in  the  6th  measurement  interval. 

Pi*(0 

ith  correlation  coefficient  from  September 

p(0 

received  digit,  waveform. 

1973  experiment. 

P«) 

calculated  received  digit  waveform. 

p(4,  j) 

sample  of  normalized  cross  correlation  function 

P<6,  0 

sample  of  calculated  digit  waveform. 

in  4th  measurement  interval. 

PK 

bit  error  probability. 

p>*(4) 

correlation  coefficient  in  6th  measurement 

4(0 

symbol  composition  coefficient  for  the  ith 

interval  from  January  1974  experiment. 

digit. 

U) 

greatest  integer  in  x. 
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flip  most  interesting  results  of  t ho  lose  nrcli  is  the* 
agreement  obt uitied  between  the  experimental; 
data  and  the  predictions  of  t his  theory. 

((’)  Because  the  physical  channel  is  distinctly 
time-varying,  the  received  symbol  waveform  must 
he  continually  measured  to  maintain  the  required 
match  between  the  filter  and  the  symbol  wave¬ 
form.  This  measurement  is  made  possible  by 
including  a  known  probe  component  in  the  trans- 
ini»siott  in  addition  to  the  unknown  information 
component.  Thus,  the  communication  systems 
described  here  perform  dual  roles — channel  meas¬ 
urement  and  information  transmission. 

(lr)  The  experimental  program  has  two  closely 
related  stages.  The  first  stage  conducted  com¬ 
munication  experiments  between  fixed  sites  to 
obtain  a  statistically  significant  quantity  of 
performance  data.  Table  1  depicts  the  history  of 
the  fixed-site  communications  experiment.  Only 
the  results  of  the  recent  experiments,  M5  and  M6, 
are  presented  here.  Previous  experiments  obtained 
similar  results,  but  tliey  were  based  on  very  few 
data.t-'1)  Sections  11,  III,  and  IV  of  this  urticlo 
detail  the  operation  and  performance  of  the  M5 
and  Mli  experiments. 

(1)  The  second  stage  of  the  experimental  program 
was  to  determine  the  applicability  of  the  fixed-site 
results  to  the  submarine  communication  problem. 
For  such  an  application  to  be  possible,  the  acoustic 
medium  must  be  stable  in  space  as  well  as  time. 
That  is,  the  probe  measurement  of  the  channel 
must  remain  valid  under  spatial  displacement 
of  the  submarine  platform.  Consequently,  n  careful 
measurement  of  the  spatial  stability  was  made,  us 
described  in  Section  V. 

((')  The  union  of  the  results  from  the  fixed-site 
communication  experiments  and  the  spatial- 
stability  measurements  solidly  indicates  the  feasi¬ 
bility  of  coherent,  matched-filter  techniques  for 
submarine  communication  systems.  'Phis  conclu¬ 
sion  and  others  are  presented  in  Section  VI.  An 
application  of  these  techniques  in  conjunction 
with  a  randomized  transmission  format  is  also 
given. 

II.  TRANSMISSION  FORMAT 

((.')  Because  of  the  dual  nature  of  the  transmission 
and  the  implementation  of  the  receiver  processing, 


Table  1.  Summary  of  fixed-site  communication 
experiments.  (Table  unclassified.) 


Experiment 

Year 

Disunite  (nm) 

Duration  (hr) 

Ml 

1968 

7 

24 

M2 

1969 

7 

24 

M3 

1969 

7 

24 

M4 

1970 

42 

4 

M5A 

1972 

7 

2214 

M5B 

1972 

7 

1500 

M6A 

1973 

42 

168 

M6B 

1973 

42 

984 

M6C 

1973 

42 

140 

the  signals  transmitted  by  the  communication 
systems  have  u  complicated  format.  One  part  of 
the  signal,  the  probe  component,  allows  measure¬ 
ment  of  the  channel;  the  other,  the  information 
component,  contains  the  information.  Subsequent 
parts  of  this  section  describe  these  components  in 
detail  und  explain  the  choice  of  transmission 
format.  Section  III  on  receiver  processing  com¬ 
pletes  the  explanation.  The  probe  and  information 
components  ure  transmitted  interleaved  in  time 
(time  multiplexed)  as  described  in  Section  Iln. 

A.  BASIC  SIGNAL  ELEMENT,  THE  DIGIT 

(U)  Both  components  of  the  transmission  are 
composed  of  a  succession  of  biphase  modulated 
elements,  ci  lied  digits.  The  simplest  example  of 
a  digit  waveform  is  a  rectangular  carrier  pulse. 
The  probe  ami  information  components  consist  of 
modulated  digits  in  a  prescribed  (or  in  the  case  of 
the  information  component,  almost  prescribed) 
order.  In  the  subsequent  discussion,  no  conflict 
will  arise  if  the  digit  waveform  is  assumed  to  be 
a  rectangular  carrier  pulse,  although  in  practice, 
some  amount  of  bundspreading  of  the  digit  is 
desirable,  as  suggested  below. 

(U)  For  a  rectangular  carrier  pulse,  the  pulse 
duration  and  bandwidth  are  inversely  related; 
that  is,  once  the  duration  is  specified,  the  band¬ 
width  is  also  fixed.  This  constraint  is  undesirable 
in  practice,  since  the  choice  of  digit  duration  and 
digit  bandwidth  is  based  on  different  con¬ 
siderations.  Typically,  the  digit  duration  is 
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Figure  1m.  Digit  waveforms:  rectangular  elicit.  (Figure 
uncln^ilicd.) 

«  win 


•T/J 


Figure  II).  Digit  \vm\ cfnnus :  haudsprcad  digit  with  three 
silliclcmcnt,*.  (Figure  unclassified.) 

selected  on  the  basis  of  nitihijmtli  conditions, 
while  tin1  digit  bandwidth  is  based  on  the  trans¬ 
duce!'  bandwidth. 

((’)  In  gmtornl,  the  problem  is  to  increase  the  elicit 
bandwidth  beyond  that  of  a  rectangular  digit  of 
t*(|iml  duration.  Figure  1  shows  the  results  of 
bandspreiidiug  of  the  tli^ri t  waveform  by  compari¬ 
son  with  a  rectangular  pulse  and  its  spectrum  in 
Fig.  la.  In  each  [tart  of  Fig.  1,  i/(t)  is  the  digit 
waveform  and  £/(/)  is  the  digit  spectrum.  Note 
that  as  the  number  of  subelements  in  the  digit 
waveform  is  increased,  the  bandwidth  increases 
proportionately.  The  digit  waveform  in  Fig.  lb 
was  used  in  the  Ml,  M2,  M-i,  and  M4  experiments, 
and  the  waveform  in  Fig.  1c  was  used  in  the  MOB 
and  MGC  experiments. 

H.  PIlOtlE  COMPON’EN'T 

((')  The  probe  component  of  the  transmission 
allows  the  receiver  to  measure  the  channel  digit 


um 


•r/j 


Figure  lc.  Digit  waveforms:  bandspread  digit  with  seven 
subelcmcnts.  (Figure  midassilicd.) 

response.  Because  the  information  symbols  are 
composed  of  combinations  of  digits,  a  filter 
matched  to  the  symbols  can  he  formed  from  the 
digit  response.  The  probe  component  must  be 
constant  and  known  to  the  receiver  if  this  measure¬ 
ment  is  to  be  successful.  This  is  in  contrast  to  the 
information  component,  described  later  in  Section 
lie,  which  is  inherently  variable  and  unknown. 

(U)  One  method  of  measuring  the  channel  digit 
response  is  to  transmit  a  single,  isolated  digit  us 
the  probe  component,  as  shown  in  Fig.  2u.  Vacant 
time  intervals  must  he  included  before  and  after 
the  probe  digit  to  allow  for  the  time  smear  of  the 
acoustic  channel.  The  disadvantage  to  this  probe 
structure  is  that  energy  is  transmitted  for  only  a 
small  fraction  of  the  probe  duration.  Consequently, 
the  receiver’s  measurement  of  the  channel  is  noisy. 

(U)  A  belter  choice  for  the  probe  component  is  to 
form  it  from  an  integral  number  of  periods  of  a 
linear  maximal  sequence.  This  allows  more  energy 
to  be  transmitted  in  the  probe,  with  a  correspond¬ 
ing  improvement  in  the  measurement  of  the 
channel  digit  response. 

(U)  Let  u(t)  be  the  transmitted  digit  waveform  of 
duration  T.  Let  b(i),  equal  to  cither  +1  or  —  1 
[6( /’)* |  ;fc  1 1  ],  he  the  values  of  a  Km  long,  linear, 
maximal  sequence.  The  orohe  component,  mv{l), 
of  the  transmission  i-  obtained  by  modulating 
successive  digits  with  KlM  periods  of  the  sequence: 

nip(l)  2  b(i)u(t  —  iT)  (1) 

where 

i\\m  =  K^K,m  (2) 
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The  transmitted  probe  component  will  be  time- 
limited  to  n  duration  Ttm—Xlm  T,  because  of  the 
limited  duration  of  the  individual  digits.  A  probe 
component  of  (his  form  is  shown  in  Fig.  2b  for 
K,m—  1  and  /v*,=  15. 

(U)  Because  of  the  assumed  linearity  of  the 
acoustic  channel,  the  received  probe  xr(1)  is  a 
superposition  of  channel  digit  responses,  p(t). 

.Vfm-l 

/h(0=  Z  bd)p(t-iT)  (3) 

i*»Q 

In  general,  time  smearing  of  the  channel  will  make 
the  duration  of  the  received  probe  longer  than  TtUi. 

(U)  The  primary  objective  of  the  probe  component 
is  to  allow  measurement  of  p(t) ;  and  consequently 
/,.(/),  as  indicated  in  Eq.  (3),  must  be  processed 
further  to  be  useful .  A  simple  digital  filter  can  be 
applied  to  sv(t)  to  obtain  u  waveform  p(t),  which 
approximates  the  channel  digit  response  p(t).  To 
understand  the  derivation  of  ,ne  approximate 
channel  response  /;(/),  consider  the  following 
analysis  for  an  infinitely  long  periodic  probe 
component,  /,."(/) 

J’p" ( t)  —  2  b(i)p(t—iT)  (4) 

la*® 

Define  p”U),  the  calculated  channel  response,  as 
follows: 

]>*«)  -  g''  6-°2—  Mt+jT)  (5) 

Note  that,  although  /,."(f)  is  of  infinite  extent, 
the  summation  for  p*(0  involves  only  Ntm 
terms. 

(U)  If  Eq.  (4)  is  inserted  into  Eq.  (5),  the  following 
expression  is  obtained: 

A'ln,-1  h(  i)~ H 

/'”(«)=  z  z  b(i)p[t-a-j)T) 

(6) 

After  a  change  of  variables,  k=i—j,  Eq.  (6) 
becomes 

Z  P(t-kT)  2  b(j+k) 

;  -o  £ 

(7) 


Figure  2.  Probe,  component  structures.  (Figure  Rnd  figure 
caption  each  classified  Confidential.) 


Because  of  the  properties  of  linear  maximal 
sequences,  the  summation  over  j  is  known: 

"g’  b{j+k)A£±  •  K,m  (8) 

for  k—0,  K„,  2 Km,,  .  .  . 

Ng‘  ^±1  b(j+k)=0  (9) 

otherwise.  With  the  above  results  substituted  into 
Eq.  (7),  the  following  simplified  expression  1- 
obtained : 

£•(,)-_  *!a±l  Kfm  £  V(t-Kn  •  t T)  (10) 

That  is,  p“(t)  is  a  periodic  (with  period  7’*,= 
Kk-T)  replica  of  p(t)  with  a  processing  gain  of 
#tm-(A\y+ 1)/2.  Consequently,  with  an  infinitely 
long  probe  component,  Eq.  (5)  allows  the  exact 
measurement  of  the  chunnel  digit  response. 
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l(’)  Thi‘  probe  component,  s,.(t),  for  tin  actual 
commmiicaliou  system  cannot  have  infinite  ex¬ 
tent,  us  does  However,  a  periodic  exten¬ 

sion,  r,.Eit),  can  be  constructed  from  fv(t)  to 
yield  an  approximation  to  p(t).  Assume  that  most 
of  the  energy  of  sv(t )  is  contained  in  the  time 
interval  (/„,  t„+  T,m).  Consider  the  infinite  periodic 
extension,  rvK(t),  of  xv(t)  in  this  interval: 

.'pe(0  =  /,.(/) ;  ta<t<t0+T,m  (11) 

—  xv(t  —  k'Ti m) ;  /0+£  T,u>  fs  (&+ 1)  T,m 

(12) 

If  A’Im>l  and  the  time  duration  of  />(/)  is  small 
relative  to  the  sequence  durution,  T*,—  A'*,  T,  then 

/pK(0=/e”(0  03) 

So  an  approximation  lo  pit)  is  calculated  as: 

•Vim- 1  hi ; )  _u  i 

/»(/>  jg  U2 r  /eE(0;  t0<t<Tn  (14) 

■  K,mp(ty,  t0<t<t0+Tt!l 

(.5) 

That  is,  p(t)  is  approximately  equal  to  pit)  multi¬ 
plied  by  a  constant.  '1'his  constant  gives  a  signifi¬ 
cant  processing  gain  against  noise. 

((')  In  the  M5  and  M6  systems,  Km  was  selected 
during  the  system  design  stage  lo  be  15,  so  Tm 
exceeded  the  expected  smear  of  pit)  by  u  factor  of 
at  least  three.  The  constant  Ktw  was  selected  to 
be  four  so  as  to  support  the  approximation  indi¬ 
cated  by  Kq.  (13).  With  these  parameter  values, 
the  use  of  the  coded  probe  component  yielded  a 
9-db  processing  guin  over  a  single  digit  probe 
repeated  A’fm  times.  No  adverse  effects  due  to  the 
approximation  were  observed. 

C.  INFOHMATION  COMPONENT 

(C)  The  information  component  can  be  formed  in 
one  of  three  ways:  it)  biphase  modulation  of 
single  transmitted  digits;  (2)  biphase  modulation 
of  groups  of  digits,  where  the  intragroup  structure 
is  constant:  or  (.?)  variable  symbol  modulation, 
in  which  biphase  modulation  of  groups  of  digits, 
where  the  intragroup  structure  varies  from  one 
symbol  to  the  next,  is  carried  out. 


((')  In  the  subsequent  discussion,  the  difference 
between  a  digit — described  previously — and  a 
symbol  is  important.  A  symbol  is  a  digit  or  group 
of  digits  used  to  carry  a  single  bit  of  information. 
The  objective  of  the  receiver  is  to  determine 
symbol  values,  not  digit  values;  thus,  the  difference 
between  the  two  methods  of  forming  the  informa¬ 
tion  component  lies  in  the  structure  of  the  symbol 
waveform. 

1.  Single  Digit  S'l/mbols 

(U)  The  simplest  construction  of  the  information 
component  assigns  a  single  digit  to  each  symbol, 
A  digit  transmitted  at  0  phase  represents  a  binury 
one;  a  digit  transmitted  at  180-deg  phase  repre¬ 
sents  a  binary-minus  one.  This  technique  wus 
used  in  the  Ml,  M2,  M3,  and  M4  systems. 

(C)  In  all  of  the  systems  considered,  the  number  of 
digits  in  the  information  component  was  made 
equal  to  the  number  of  digits  in  the  probe  com¬ 
ponent.  When  the  probe  component  contained 
I\ io,  repetitions  of  the  K digit  linear  maximal 
sequence,  the  information  component  was  com¬ 
posed  of  A',m  repetitions  of  K „  information- 
bearing  digits.  This  symmetric  strut  ture  is 
arbitrary,  i.od  it  effectively  sets  the  information 
and  probe  component  energies  equal.  In  general, 
the  structures  of  the  probe  and  information  com¬ 
ponents  are  independent,  subject  only  to  the 
constraint  that  both  be  composed  of  the  same  digit 
waveform  Thus,  the  A*(!i  transmitted  information 
component  m,(k,  t)  is  given  bv 

nti{k,  t)  2  X  c(A'„g-  k+j)n[t— (A’ia-  i+j)T] 

i»(i  i~o 

(10) 

where  c(/i)tj  ±  1  J  is  the  at  It  information  bit. 

(lT)  Figure  3a  depicts  an  information  component 
using  a  single  digit  per  symbol.  The  diagram  shows 
A’„,  symbols  corresponding  to  the  information 
stream  -1  +  1 -1  +  1  +  1,  -1-1  +  1 -1-1, 

—  1  + 1  + 1  + 1  +  I .  Only  one  repetition  ( K,m  -  1 )  is 
shown  for  simplicity.  In  actual  practice  A",  „  is 
taken  as  being  larger  than  one  to  increase  the 
energy  content  of  each  symbol  and  to  provide 
limited  time  diversity. 

2.  Multiple  Digit  Symbol * 

((’)  If  each  information  symbol  is  composed  of 
only  one  digit,  as  described  above,  then  potential 
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tradeoff  opportunities  for  system  design  are 
eliminated.  For  example,  the  number  of  digits  in 
the  sequence,  Km,  and  the  digit  duration,  T,  must 
be  selected  so  that  the  time  duration  of  the 
received  digit,  pit),  is  less  than  the  sequence 
period,  7’^,  in  order  to  properly  measure  p(t). 
Because  of  memory  requirements  in  the  receiver, 
the  sequence  period,  Tn,  should  b<>  as  short  as 
possible.  However,  the  resulting  digit  duration,  T, 
nun  be  too  small  to  provide  adequate  symbol 
energy  and  to  avoid  intersymbol  interference. 

t( ')  One  solution  to  the  above  problem  is  to  simply 
allow  more  than  one  digit  per  symbol,  keeping  the 
digit  duration,  T,  and  Km  constant.  In  the  M5 
and  Mt>  systems,  three  digits  were  used  for  each 
information  symbol,  giving  a  threefold  increase  in 
both  symbol  energy  and  duration,  with  no  increase 
in  sequence  duration,  TK. 

((')  The  composition  of  the  information  symbol  is 
arbitrary,  as  long  as  it  is  composed  of  linear 
combinations  of  digit  waveforms.  Let  y(t)  be  a 
symbol  that  is  taken  us  being  composed  of  Km 
digit  waveforms, 


c/(f)-  qii)n(t-iT) 


Figure  3.  Information  component  structures.  (Figure  and 
(17)  figure  caption  each  classified  Confidential.) 


where  f/(  / )« !  ±  1 1 ;  /  =0,  .  .  .  ,  K* m  —  1 .  The  choice 
of  digit  reparations  of  multiples  of  T  sec  and  unity 
magnitudes  for  the  c/(»)’s  is  not  essential  in  general, 
but  it  was  made  for  ease  of  generation  and  pro¬ 
cessing  of  the  signal.  If  pit)  is  the  received  digit 
response,  then  a  filter,  hit) ,  matched  to  the 
received  symbol  is 

h(t)  -  ^  qd)p(t-iT)  (18) 

by  the  assumed  linearity  of  the  channel.  Thus,  u 
filter  matched  to  a  specified  symbol  can  easily  be 
constructed  from  the  received  digit  waveform  p(t) 
and  the  symbol  composition  given  by  </(0),  .  .  .  , 
'/(A'„m-  1). 

(U)  Figure  3b  depicts  un  information  component 
with  A'„„  =  3  digits  per  symbol  and  with  2o=L 
Vi  —  I ,  and  q<=  —  l.  The  information  values  shown 
are  —  1  +  I  —  1  + 1  -f- 1 .  Again,  A,m=  1  for  the  sake 
of  simplicity. 

CONFIDENTIAL 


S.  Variable  Symbol  Modulation 
(C)  If  more  thun  one  digit  is  contained  in  each 
symbol,  us  described  above,  the  opportunity  to 
alter  the  symbol  composition  from  one  time 
interval  to  the  next  arises.  For  example,  consider 
K„  consecutive  symbols  having  Km  digits  per 
symbol.  The  jth  symbol  waveform  is  specified  b\ 
K,m  coefficients,  y(0)  .  .  .  q(Km  —  1),  q(i)t{±\}, 
giving  a  possibility  of  2K,m  different  symbol 
waveforms.  There  is  no  requirement  that  succes¬ 
sive  symbols  have  the  same  composition,  so  one 
could  define  K,t  symbol  waveforms  y(j,t)  j—0  .  .  . 
K„- 1,  each  with  different  composition, 

Kmm  —  i 

t i(j,  0=  ^  vO.  j=0  .  .  .  Ktt-\ 

(19) 

to  be  transmitted  Tm=K,mT  see  apart.  Because 
of  the  assumed  linearity  of  the  channel,  a  filter 
matched  to  each  symbol,  h(j,t),  ean  easily  be 
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constructed  from  pit).  The  *th  transmitted 
probe  eoinponent  m,(L\t)  is  then 

r>n(k,  t)  ^2  2  tlKu  ■  k+j) 

t-n  j^n 

X  g(j,  t-j  .  Tm-iT„)  (20) 

The  teehnique  in  which  the  symbol  composition 
is  permitted  to  vary  from  one  symbol  to  the  next 
is  called  variable  symbol  modulation. 

(C)  Variable  symbol  modulation  was  incorporated 
into  the  M5  and  M6  communication  systems  as  a 
means  of  facilitating  the  experimental  measure¬ 
ment  of  system  performance.  In  Section  VI h, 
variable  symbol  modulation  is  shown  to  lx*  an 
important  component  of  a  detection-resistant 
transmission  format.  In  the  M5  and  M6  systems, 
the  product  of  K,„,  and  Klt  was  selected  to  equal 
A«q,  specifically,  A’o— 5.  Under  these 

conditions,  the  probe  and  information  components 
can  be  made  identical  by  a  proper  choice  of 
qii.j)  and  the  information  values  e(i). 

(C)  Let 

'/O'.  i)=b(K,m  •  j+i); 

*=°  .  .  .  K.m—  1  and  j= 0  .  .  .  Kta—  1  (21) 

where  6(f),  >’= 0  .  .  .  K „ —  1  are  the  biphase 
sequence  coefficients  described  in  Section  IIb. 
Then  if  r(j)  —  1 ,  for  all  j,  the  transmitted  informa¬ 
tion  component  is 

»«•(*.  o  =25'  s'  e(K„  ■  k+j)  q(j,  „) 

/“  0  w=0 

X  i<[t—(iT„+jTm+nT)]  (22) 
“  JC  $2  b(j)ii[t-(K,Q  -  i+j)T\  (23) 

S  b(j)"(t—jT)  (24) 

"  hich  is  exactly  the  same  as  the  prol>e  component 
given  by  Eq.  (3). 

(U)  The  advantage  of  the  above  choice  of  Km, 
and  e{i)  is  that,  since  the  probe  and 
infonnation  components  are  identical,  a  periodic 
transmission  is  a  valid  system  transmission.  In 
measurement  of  the  communication  performance 
of  the  .\I5  and  Mt>  systems,  such  a  periodic  wave¬ 
form  was  transmitted.  ’I  hLs  simplified  the  structure 
of  the  transmitter  and  eliminated  the  need  for 


synchronization  of  the  receiver  during  the  lengthy 
evaluation  process. 

(U)  Figure  3c  depicts  an  information  component 
with  variuble  symbol  modulation.  The  informa¬ 
tion  values  being  transmitted  are  —I,  +1,  —  l, 
+ 1 ,  + 1 ,  and  Ktm  =  1 . 

D.  MULTIPLEXING  OF  PROBE  AND  INFORMATION 
COMPONENTS 

(C)  The  transmission  in  each  of  the  matched- 
filter  communication  systems  consists  of  a 
multiplex  of  the  probe  and  information  com¬ 
ponents.  In  choosing  the  multiplexing  technique, 
the  particular  purpose  of  the  probe  component, 
must  be  considered.  For  example,  frequency 
multiplexing  would  be  inappropriate  in  the  acous¬ 
tic  medium,  because  the  digit  response  measured 
in  one  frequency  bund  would  not  be  valid  in 
another  frequency  band.  Two  suitable  multi¬ 
plexing  techniques  are  time  multiplexing  and 
phase  multiplexing.  The  M5  and  MG  systems  used 
time  multiplexing,  as  described  below.  Later,  in 
Section  VI b,  an  application  of  phase  multiplexing 
is  presented. 

(C)  In  the  M5  and  M6  systems,  the  transmission 
m(0  ,vus  composed  of  a  simple  time  multiplex  of 
probe  und  information  components: 

-f  <o 

»(<)=  X!  mP(t-2Tla  ■  O+Mili,  /-(2H-l)7\m] 

Is-® 

(25) 

where  T,m  is  the  probe  (and  information)  compo¬ 
nent  duration  and  Ttm=Ktni-Kmi-T. Thus,  the  trans¬ 
mission  consists  of  alternating  probe  and  info:- 
mation  components,  with  each  successive  in- 
formation  component  representing  different 
information  bits.  Figure  4  depicts  a  typical  trans¬ 
mission  from  the  M5  and  Mo  systems.  The 
information  values  for  the  information  component 
shown  in  expanded  form  are  — 1,  -f  1,  — l,  -f-i, 

III.  RECEIVER  PROCESSING 

(C)  The  receiver  in  the  M5  and  MG  communi¬ 
cation  systems  performs  a  duul  role.  The  first 
part  of  this  dual  role  is  to  measure  the  existing 
channel  digit  response  from  the  probe  component 
of  the  transmission.  Several  additional  measure¬ 
ments  are  made  from  the  probe  in  order  to  evul- 
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milt'  system  performance.  The  second  part  of  the 
receiver's  role  is  to  determine  the  transmitted 
symbol  values.  In  general,  time  synchronization 
is  required  to  direct  the  probe  component  of  the 
transmission  to  the  probe  processing  section  of 
the  receiver,  etc.  This  synchronization  is  discussed 
ip  Section  1 1  Id,  while  the  following  discussion 
(particularly  Sections  Min  and  IIIc)  assumes 
that  correct  synchronization  has  already  been 
obtained. 


A.  PRELIMINARY  PROCESSING 

(lT)  Signals  from  the  receiver  hydrophone  are 
transmitted  through  linear,  fixed-gain  amplifiers 
to  the  processor,  the  gains  of  these  amplifiers 
being  selected  so  that  no  clipping  occurs  during 
normal  operation.  At  the  processor  site,  a  fixed 
analog  filter  centered  on  the  carrier  wave,  420  hz, 
and  with  nominal  bandwidth  of  100  hz  removes 
much  of  the  ambient  noise  from  the  signal.  The 
output  of  the  filter  is  amplified  by  another 
fixed-gain  amplifier  and  applied  to  the  digital 
processor  nnulog-to-digital  converter. 

(U)  The  analog  converter  in  the  processor  handles 
signals  in  the  ±  I -volt  range  and  bus  either  8 
bit  (Mo)  or  9  bit  (MCA,  MOB,  M6C)  plus  sign 
resolution,  with  the  equivalent  dynamic  ranges 
being  48  and  54  ill),  respectively.  System  gains 
are  adjusted  so  (hat  normal  operation  occurs  ut 
roughly  10  dl>  below  the  converter  limit.  Local 
shipping  noise  does  cause  overloads  to  occur 
throughout  the  system,  but  these  overloads  are 
infrequent  and  are  eusily  recognized. 

(LT)  The  sampling  of  the  bandpass  input  signal  is 
controlled  by  a  high-stability  oscillator  (at  least 
I  part  in  !();)  set  for  h  sampling  rate  /c(,  of  four 
times  (1080  hz)  the  currier  frequency.  The 
samples,  are  digitally  processed  to  obtain 

complex  low-pass  samples,  /(»),  equivalent  to  the 
original  bandpass  waveform. 


*(»)-■  J-  £  (— \)‘s'(2Kt9  •  i+2j)  (i  even) 

Acl,  JTo 

(2«) 

'(>■>  K  £  (-1)' 

A.n, 

X  /Wcp(i— l)+2;+lJ  (iodd)  (27) 


In  the  above  equations,  Kct  is  the  compression 
factor.  The  samples  /{>)  with  even  indices  can  be 


i - - «»  - .  i 


Figure  4.  Typical  trim-mission  from  the  M.i  and  Mfl 
systems.  (Figure  classified  Confidential.) 


associated  with  the  real  component  of  the  signal, 
und  the  odd  indices  can  be  associated  with  the 
imaginary  component.  The  equivalent  low-pass 
bandwidth  IT,*  of  either  component  is  given  by 


W 


■p 


ftlk_ 


(28) 


where  /ci|,  is  the  sampling  frequency. 

(U)  Two  additional  constants  will  be  helpful  in 
working  with  the  compressed  samples  i\i)  of  the 
input  waveform.  The  first  is  the  number  of  samples 
per  digit,  N^,  given  by 

A^%—  (29) 

IV  op 


The  sec  ond  constant  is  the  number  of  compressed 
samples  per  sequence  period  Nm: 

(30) 


As  will  be  seen  shortly,  is  an  important 
constant  in  determining  the  memory  requirements 
of  the  receiver. 


B.  PROBE  PROCESSING 

(C)  The  primary  objective  of  the  probe  processing 
is  to  determine  the  channel  digit  response,  p(t). 
A  secondary,  but  important,  objec  tive  of  the  probe 
processing  is  to  measure  the  basic  transmission 
characteristics  of  the  medium  to  aid  in  evaluation 
and  understanding  of  the  system  performance. 
All  of  the  communications  systems  included 
measurement  of  the  widebund  signal  power  and 
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noise  power,  SP  and  NP,  while  several  specialized 
measures  of  the  time  smear  mid  time  stability  of 
the  channel  are  also  derived  from  p(t). 

(C)  The  probe  processing  operates  on  the  probe 
component  of  the  transmission.  At  the  input 
sigual-to-noise  ratios  below  0  db  for  which  the 
systems  were  designed,  the  energy  contained  in  a 
single  probe  component  is  inadequate  for  obtaining 
a  satisfactory  estimate  of  the  channel  digit 
response,  /;(/).  Because  of  this,  coherent  integra¬ 
tion  of  Km i  probe  components  is  performed,  with 
the  time  interval  over  which  it  is  performed  being 
called  the  measurement  interval.  To  reduce  the 
receiver  memory  requirements,  this  integration  is 
performed  in  two  steps. 

(C)  The  first  step  is  to  integrate  over  the  Km 
sequence  periods  contained  in  a  single  probe  com¬ 
ponent  to  yield  a  single  vector  of  length  Nm.  Let 
aAk,  i)  be  the  /th  element  of  this  vector,  derived 
from  the  4th  probe  component.  That  is, 

«p(*.  *)—  S  42  k-  N „  •  K,m+j  ■  N„+i) 

0<k<Kal-,  0 <i<Nmt  (31) 

The  above  equntion  assumes  that  the  receiver  has 
been  synchronized  so  that  j-(0)  .  .  .  r(NK  K,m 
—  1)  are  samples  of  a  probe  component. 

((')  The  second  step  in  the  integration  of  the  Kml 
probe  components  is  to  derive  the  signal  vector 

aAi): 

Op (*)—  aAk,  i)  (32) 

The  signal  vector  aH(»)  is  the  coherent  integration 
of  K,m  Kn,  sequence  periods  and  is  the  starting 
point  in  the  determination  of  the  channel  digit 
response,  At  the  same  time,  another  vector 
of  points  is  also  calculated. 

n(i)=  ^  (— l)*oP(t,  t)  (33) 

The  noise  vector,  «(<),  eliminates  the  signal  con¬ 
tribution  by  alternating  the  sign  in  its  summation, 
and  it  is  the  basis  for  the  noise  power  measure¬ 
ment. 

1.  Channel  Digit  Response  Measurements 
(U)  As  mentioned  above,  the  primary  objective 
of  the  probe  processing  is  to  derive  the  approxi¬ 


mate  channel  digit  response  p(t).  Let  /,(*,/)  be 
the  sampled  data  representation  of  p(t)  obtained 
in  the  irth  measurement  interval.  From  a  discrete 
version  of  the  derivation  of  Section  IIb,  the 
following  result  is  obtained : 

P(*.  <)=*§'  -(I)2~  Mi+j  ■  Ard»);  0<i<-V., 

(34) 

where  the  index  is  taken  modulo  Nm.  Tho 

index  k  on  p(k,  i)  is  introduced  ut  this  point  to 
facilitate  subsequent  description;  it  indicates  that 
p(k,  /')  is  derived  from  the  4th  measurement 
interval  of  km,  probe  components. 

(U)  An  oscilloscope  display  of  the  p(k,  i)  vector  is 
provided,  and  it  is  updated  during  each  measure¬ 
ment  interval  by  the  processor.  Unfortunately,  the 
Nn  points  in  the  display  are  too  many  to  be 
retained  during  a  long-term  experiment;  however, 
two  useful  measures  of  the  multipath  situation  are 
recorded  for  evaluation  of  system  performance. 

(U)  The  first  multipath  measure  is  the  autocor¬ 
relation  of  the  channel  digit  response,  p(k,  i),  taken 
at  multiples  of  the  digit  length,  Nai. 


K;-,(k,  P(k,  j) 

X  P(k,  j+iN«)i  i=0  ...  6  (35) 

where  the  index  j  +  /Ara,  is  taken  modulo  Af„. 
These  measurements  are  indicative  of  the  time 
smear  of  the  channel,  and  they  are  important  in 
measuring  the  extent  of  intersymbol  interference. 
For  example,  if  the  channel  were  idcul  (broad 
band)  then  p(k,j)  and  j>(k,j+iNiK)  would  have  no 
overlap  for  <  >0.  Then  //- -(A*,  iNae)  would  be  zero 
for  />0.  Thus,  as  //"(A-,  nVds)  with  /> 0  becomes 
large  relative  to  0),  the  time  smear  of  the 

channel  increases. 

(U)  The  second  multipath  measure  is  the  cross 
correlation  between  two  consecutive  measure¬ 
ments  of  fi(k,t).  Let  p{k—  l,f),  p(k,t)  be  two  digit 
response  measurements,  based  on  consecutive  but 
nonoverlapping  measurement  intervals.  Then,  the 
normalized  cross  correlation  coefficient,  p(k),  be¬ 


tween  the  two  measurements  at  zero  offset  is  given 

by 


2b  i>(k—  1,  i)p(k,  i) 

P(*)=-"0 


V*;;(fc—  l,  0)Rn(k,  o) 


(36) 
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Tin*  mufriiit iido  of  p(k)  is  always  less  than  one.  If 
(In'  nroii't  ir  channel  «i'iv  I  nil  y  l  into  invariant 
(and  thi'ii'  were no  iium*),  (lien  p(h  would  equal  one. 
If  the  eluinnel  wa>  complete!)  unstable  so  that 
the  two  parameter*  i )  and  jt(k,  <)  had 

not hin*;  in  rninuion,  then  pile)  would  have  a  mean 
Value  of  zero. 

Mnixiiinnml  of  Slijnal  anil  A’o/w  Power* 

( l')  The  wideband  signal  and  noise  powers.  SP  and 
Nl*.  are  measured  to  allow  evaluation  of  system 
performunee.  Let  *(/')  he  the  signal  eontribution  to 
the  dll  'ample  in  a  single  received  sequence  period, 
and  lei  «(»)  he  the  noise  eontribution  to  the  dll 
sample  in  a  single  received  sequence  period.  The 
noise  contributions  n(i)  are  assumed  to  be  nor¬ 
mally  di'tributed  with  zero  mean  and  variance 
<rN-‘ and  are  also  assumed  to  be  independent.  Then, 
if  the  channel  is  stable, 

cf,.(i)-A'(111  •  K„,i*0')+  2  n(j  •  N«,+i) 

(37) 

Then  av(i)  w  ill  be  of  mean  KlmKm,sfi)  and  variance 
K,m-Kni,as-.  The  average  signal-to-noise  ratio 
(mean-squared  to  variance  ratio)  on  a, .(/),  SNR„UI 
is  then 


SNR m,=K„nKml  [s'(0/  (38) 

■SO  Oh 

The  constants  Ks„  and  Km  are  selected  at  the 
time  of  system  design  so  that  the  noise  contribu¬ 
tion  to  «,.(<’)  can  be  neglected.  The  approximate 
signal  power,  SP,  is  given  by 

SP--=A£'  MOI2  (39) 

i-0 

^(Ktm-K^y^  Hi)]1  (40) 

i=0 

which  is  the  signal  energy  over  Kmt  probe  com¬ 
ponents. 

((.’)  If  the  medium  is  stable  over  the  measurement 
interval,  the  signal  eontribution  to  the  n(i) 
vector  is  zero.  Thus,  nli)  has  a  zero  mean  and  u 
variance  K,w-Kwi-os--  Let 

NF  -  A£'[«(i)]J 

i  *0 
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Then  NP  is  x*  distributed  with  Nm  degrees  of 
freedom.  In  the  systems  discussed  here,  jV„>>50 
and  the  distribution  on  NPcun  be  approximated 
by  a  (Jaiissiun  distribution  with  meun 
Abc  os*  and  vurianee  ’l'he  mean- 

squared  to  variance  rutio  on  NP  is  high,  so  NP 
is  a  good  estimate  of  the  noise  power  in  Kmt 
probe  components. 


(U)  The  input 
defined  as 


■siguul-to-noi.se  ratio,  SNRln,  is 


SNR,. 


fn'Ojdt 


(42) 


where  the  integration  is  carried  out  over  the  probe 
measurement  intervul.  From  the  sampling 
theorem, 

yv<  t)dt=^w  ai««)i'  (43) 

fnHt)dl=^Un{i)Y  (44) 


where  the  summation  extends  over  the  Nm,= 
N^K, m-Kml  samples  in  the  probe  measurement 
interval.  Because  the  noise  samples  n(i)  are 
independent  and  have  zero  mean, 


E\S*M)dt)  =  ±.N„-K, (45) 
=2Tf>  (NP)  (46) 


Similarly,  the  summation  in  Eq.  (4d)  can  be 
simplified  if  the  channel  is  invariant  over  the 
measurement  interval: 

2 \y  K'“  •  K“"  [8(i)12  (47) 

1  SP 

2  IF  Ktm  ■  Kuo  ^  ; 

Substituting  Eqs.  (46)  and  (48)  into  Eq.  (42) 
gives 


Thus,  the  input  signal-to-noise  ratio  can  be  deter¬ 
mined  from  the  probe  measurements  of  SP  and 
NP.  Although  this  measurement  is  determined 
from  the  probe  portion  of  the  measurement 


ak  — i 
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interval ,  it  is  also  applicable  to  the  information 
component  by  the  assumed  invariance  of  the 
medium. 

C.  INFORMATION  COMPONENT  PROCESSING 

((’)  The  objective  of  the  information  component 
processing  is  to  make  correct  decisions  on  the 
transmitted  symbol  values.  These  decisions  can  be 
scored  against  known  answers  in  the  case  of  a 
periodic  transmission.  For  demonstration  pur¬ 
poses.  the  received  values  can  be  printed  us 
characters  on  a  Teletype.  In  either  situation, 
decisions  are  made  on  the  basis  of  thresholds  on 
the  outputs  of  a  matched  filter. 

/.  Matched-filter  Operation 

((')  The  first  step  in  the  processing  of  the  Art h 
information  component  is  to  form  a  coherent 
average  adit,  i)  of  the  information  component: 

»)—  S  1 )Ktm-Nn+jN„+i]; 

0  <i<Nn  (50) 

This  yields  a  vector  of  points  and  is  analogous 
to  the  compulation  of  the  aP(k.  i )  vector  in  the 
probe  processing.  The  «i(it,  i)  vector  contains  K,t 
information  symbols,  each  composed  of  K,m  digits 
in  a  prescribed  order.  For  the  Mo  and  M6  experi¬ 
ments.  =  :i  and  Kkt= 5,  respectively. 

((’)  To  match  filter  a, (&.'),  a  cross  correlation  is 
performed  of  ax(k,i)  with  jt(t,i),  where  />(/,/>  is  the 
digit  response  found  in  the  previous  measurement 
interval,  /.  The  decision  process  requires  samples 
of  the  matched-filter  output  taken  at  intervals  of 
one  digit  length,  ATd„  apart.  Thus,  only  cross 
correlations  are  needed  to  provide  the  required 
data  from  ut (k,i).  Specifically,  let 

£'(*»  j)  =  £  P(i+jKto)ai(k,  i);  0 <j<K„ 

1*0 

(51) 

where  the  index  /  +  jNAt  is  taken  modulo  A7*,.  The 
values  L'(k.j)  represent  a  cross  correlation  of 
<i\(k,i)  with  />(/,/)  displaced  by  multiples  of  A’d„ 
samples,  or  equivalently,  T  sec.  Note  that  the 
cross  correlation  is  performed  in  a  circular  manner, 
so  that  the  overlap  from  onp  sequence  periotl  to 
the  next  is  taken  into  account.  1 


(C)  A  natural  question  arising  from  Eq.  (51)  is 
concerned  with  how  much  synchronization  is 
required  for  the  values  L'ik.j)  to  be  taken  at 
exactly  the  correct  sampling  time.  The  answer  is 
that  nofine-qrain  synchronization  is  required.  If  the 
a,(k,i )  vector  contains  primarily  the  information 
component  (and  consequently  the  av(k,i)  vector 
contains  primarily  the  probe  component],  ns 
previously  assumed,  then  L'(k,j)  is  ulwnys  taken 
at  exactly  the  correct  time.  The  reason  for  this  is 
that,  the  transmission  provides  a  fixed  time  offset 
between  the  probe  and  information  digits.  This 
time  offset  is  completely  preserved  in  both  sections 
of  the  processor.  Thus,  if  probe  digit  A  is  separated 
from  information  digit  B  by  m  digits  in  the  trans¬ 
mission,  they  will  be  separated  by  exactly  m  {mod 
KUi)  digits  in  the  av(k,i)  und  ax(k,i)  vectors. 

(C)  The  values  L'{k,j)  are  the  response  of  a  filter 
matched  to  the  received  digit,  for  each  of  the  K„ 
digits  in  the  information  frame.  To  make  decisions 
on  the  information  symbols,  the  response  of  a 
filter  matched  to  each  of  the  symbols  must  bn 
obtained.  Since  the  information  symbols  are 
formed  from  known  combinations  of  digits,  tho 
symbol  filter  responses  can  be  found  from  combi¬ 
nations  of  the  L'(k,j). 

(C)  Specifically,  suppose  that  the  yth  symbol  in 
the  information  component  is  as  given  in  Eq.  (19). 
Then  output  of  the  matched  filter  for  the  jth 
symbol  in  the  JMt  frame,  L{k,j),  is  given  by 

/>(*, i)-As‘  '/O',  i)  s'  V\l,  (i  •  Km 

i*0  to**  0 

+£)Art«+mla,(fc,  m)  (52) 

where  the  index  (j-K,m4-/')A7d,+m  is  again  taken 
modulo  A'„.  The  summation  over  m  in  the  above 
equation  is  simply  L'(k,j-Km+i),  so  the  following 
result  is  obtained: 

L{k,j)=  S  q{j,i)L'(k,  j  ■  Km+i)  (53) 

Because  the  above  equation  is  easily,  implemented, 
the  use  of  variable  symbol  modulation  places 
very  little  additional  burden  ou  the  receiver 
processing. 

2.  Decision  Process 

(C)  To  determine  the  symbol  values,  a  zero 
threshold  is  applied  to  the  symbol  matched-filter 
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outputs  L(k,i).  If  the  filter  output  L(k,i)  is 
•neater  than  zero,  the  decision  <l(kK,,+  i)  on 
the  /tit  symbol  in  the  I-t  It  informal  ion  component 
is  I.  Otherwise  the  decision  is  —1.  In  the  MS  and 
M6  systems  K„=5,  so  each  information  com¬ 
ponent  yielded  five  bit  decisions.  These  bit 
decisions  were  printed  as  a  single  five-level 
Teletype  character  during  text  transmissions. 

(C)  More  elaborate  decision  techiques  can  be 
applied  to  reduce  the  effects  of  intersymbol 
interference.  In  the  N15  and  MG  systems,  the 
choice  of  symbol  durations  was  such  that  inter¬ 
symbol  interference  was  negligible,  and  these 
lcelntit|iies  were  not  required. 

3.  Scoring  of  Receiver  Decision* 

(U)  To  evaluate  the  system  error  performance, 
the  receiv  er  decisions  il( i)  are  scored  against  known 
correct  values.  This  can  be  (lone  as  long  as  the 
symbol  values  in  the  information  component  are 
known,  as  they  are  in  the  case  of  a  periodic  trans¬ 
mission.  The  error  count  on  the  information 
components  for  the  periodic  transmission  is  desig¬ 
nated  E,  and  is  one  of  the  parameters  computed 
for  every  Kml  probe  components. 

((’)  If  information,  unknown  to  the  receiver,  is 
being  transmitted,  it  is  of  course  not  possible  to 
score  the  information  component.  A  completely 
adequate  measure  of  system  performance  can  be 
obtained,  however,  by  performing  the  decision 
process  on  the  vector  ar(k,i)  in  a  manner  identical 
to  that  performed  on  at(k,i).  Since  the  digits 
in  the  probe  have  a  fixed  interpretation  in  terms 
of  symbols,  scoring  of  the  probe  component  is 
always  possible,  since  it  is  independent  of  the 
information  being  transmitted.  The  error  count 
ott  the  probe  component  is  designated  Er. 

((.')  When  a  periodic  transmission  is  sent,  the 
expected  values  of  the  two  error  counts  Et  and 
EP  should  be  the  same,  as  the  ae(k,i)  and  at(k,i) 
vectors  are  statistically  the  same.  On  the  other 
hund,  when  text  is  transmitted,  Et  depends  on  the 
particular  text  transmitted,  and  it  is  not  a  meas¬ 
ure  of  system  error  performance. 

D.  SYNCHRONIZATION 

((')  In  the  preceding  discussion,  synchronization 
between  the  receiver  and  received  signal  has  been 
assumed.  When  a  periodic  transmission  is  sent, 
the  probe  and  information  components  are 
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identical  and  no  synchronization  is  necessary. 
On  the  other  hand,  when  text  is  transmitted, 
synchronization  is  required  to  ensure  that  each 
component  is  processed  in  the  correct  manner. 

(C)  For  a  practical  communication  system,  syn¬ 
chronization  must  be  |>crformcd  automatically 
and  checked  continuously.  The  systems  con¬ 
sidered  here,  however,  were,  intended  primarily  to 
determine  the  feasibility  of  tin*  coherent  integra- 
tion/mulched-filter  technique.  Consequently,  very 
little  effort  was  expended  on  the  synchronization 
problem.  Indeed,  during  all  of  the  performance 
measurement  experiments  a  periodic  transmission 
was  carried  out  to  eliminate  the  problem  entirely. 
The  following  manual  synchronization  procedure 
was  used  whenever  required  for  the  actual  trans¬ 
mission  of  text. 

(C)  Synchronization  Ls  achieved  by  interaction 
between  an  operator  observing  the  oscilloscope 
display  of  p(k,  i)  and  the  processor.  The  operator 
can  delay  the  processor  time  base  for  an  integer 
number  of  digit  durations  by  typing  an  appro¬ 
priate  command.  Two  steps  are  required  to  achieve 
synchronization  with  different  transmission  being 
required  for  each  step. 

(C)  The  first  step  requires  the  basic  periodic 
transmission.  The  operator  olvserves  the  display 
of  ji(k,  1 )  and  causes  sufficient  delays  to  locate  the 
digit  response  in  the  center  of  the  display.  This 
locates  the  junction  between  the  probe  and  in¬ 
formation  components  to  be  within  an  integer 
number  of  sequence  periods. 

(C)  During  the  second  step,  a  transmission  with  a 
constant  carrier  in  each  information  component  is 
used.  If  the  receiver  is  properly  synchronized,  the 
channel  digit  response  is  a  clean  pulse,  with  zero 
levels  to  cither  side.  If  the  receiver  is  not  properly 
synchronized,  then  the  Or(k,i)  vectors  contain  the 
carrier  ns  well  as  the  sequence,  and  the  display  of 
the  digit  response  has  a  distinctive  nonzero  com¬ 
ponent.  The  operator  delays  the  processor  by 
multiples  of  Tn  until  the  proper  display  is 
obtained. 

(C)  The  preceding  procedure  was  taught  to  three 
operators  who  were  unfamiliar  with  the  total  sys¬ 
tem  operation.  Euch  was  able  to  synchronize  the 
system  within  five  minutes  or  so  under  conditions 
where  the  signal-to-noisc  ratio  was  0  db.  Synchro- 
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Pigurr  5.  Typical  multipoint  plot  from  the  M5  and  M6  systems.  (Figure  unclassified.) 


ni/.ation  itt  lower  signal-to-noisc  ratios,  is  more 
difficult,  but  it  can  be  done  by  an  experienced 
operator. 

E.  IMPLEMENTATION 

(U)  The  receiver  processing  for  both  communica¬ 
tion  systems  was  performed  by  small,  general 
purpose,  digitul  computers,  frequently  called 
minicomputers.  The  M5  systems  were  imple¬ 
mented  on  a  Digital  Equipment  Corporation 
L  INC-8  computer  with  a  4,096-word  (12  bits/ 
word)  memory.  The  M6  systems  were  implemented 
on  a  DEC  PDP-8E  system  with  an  8,192  word 
(12  bits/word)  memory. 

((  )  Both  machines  contained  a  hardware  multiply 
capability.  Such  u  feature  is  required  if  the  system 
is  to  perform  the  matched  filtering  in  real  time. 
If  it  is  assumed  that  both  the  information  and 
prolie  components  are  match-filtered  to  receive 


the  information  and  score  the  probe  component, 
the  number  of  multiplications,  naalt,  required  for 
either  component  is 


n  ,  -  V  multiplies  ...  digits 

m  matched  filter  99  component 

The  time,  Tlm,  for  a  component  is 
T,m=Kta  Nm  y* 

Jt  It 

Thus,  the  multiplication  rate,  required  is 

(56) 


(54) 


(55) 


_ Enroll 

'V.mult —  rp  ~ 
1  fm 


__  K,  •  Jtlk 

K,m .  K„ 

=788  multiplications/see 


(57) 

(58) 
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for  I  v  picul  MS  and  M6  experiment  parameters. 
'Phis  rate  is  easily  maintained  with  most  modem 
minicomputers;  however,  the  multiplication  rate 
must  nevertheless  be  considered  in  system  param¬ 
eter  design. 

(LI)  The  primary  output  of  the  processor  is  digital 
magnetic  tape — s|>ecificully,  I.IXC  format  tape. 
The  probe  measurements  and  error  counts  are 
stored  in  blocks,  with  data  from  seven  meas¬ 
urement  intervals  ami  time  information  in  each 
block.  One  tape  is  filled  even  tt  or  7  days.  Every 
time  a  tape  block  is  written,  an  entry  is  typed 
into  a  Teletype  log  for  verification  purposes.  A 
secondary,  but  important,  output  is  given  on  a 
multipoint  plotter  controlled  by  the  processor. 
This  plot  gives  the  values  of  SP  and  NP  on  a 
minute-by-minute  basis.  The  probe  error  count, 
h'H,  is  also  plotted,  along  with  a  timing  track. 
Figure  5  depicts  u  typical  multipoint  plot  and  the 
associated  scale  factors. 

IV.  EXPERIMENTAL  PROGRAM 

(C)  Five  long-term  experiments  (M5A,  M5B, 
MfiA,  MttB,  ami  MW’)  were  performed  between 
fixed  sites  in  the  Straits  of  Florida  to  evaluate  the 
effectiveness  of  coherent  integrution/inutched- 
filter  techniques.  The  dates  and  durations  for  each 
of  these  experiments  are  shown  in  Table  1.  A  total 
of  4,534  hours  (190  days)  of  valid  duta  was 
obtained  and  analyzed  over  either  7-  or  42-nnti 
distances.  Section  IVa  describes  the  differences 
among  the  five  experiments.  Section  IV i*  describes 
the  experimental  results. 

A.  DIFFERENCES  AMONG  EXPERIMENTS 

(lT)  In  each  of  the  five  experiments,  the  trans¬ 
mission  was  generated  by  the  source  off  Fowey 
Rocks  Light.*  During  the  M5  experiment,  signals 
were  received  at  u  hydrophone  located  7nmi  from 
the  Fowey  Rocks  source,  while  in  the  M6  experi¬ 
ment,  signals  were  received  at  the  Bimini  hydro¬ 
phone,  42  nmi  from  the  source.  Each  experiment 
had  more  than  one  part,  as  indicuted  by  the  third 
letter  in  the  experiment  identification;  for  exam¬ 
ple,  M6C  was  the  third  part  of  the  M6  experiment. 
The  pur|K>se  of  each  part  of  both  experiments  is 
described  below. 

*(  U)  Reference  4  descrila**  the  fixed-site  ranges  in  detail. 
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Tarle  2.  Common  features  of  the  MSA,  MSP,  MG  A, 
M6B.  anil  MOC  experiments.  ( Table  classified 
Confidential.) 

Center  Frequency.  fn  (Hz)  420 

Nominal  Bandwidth.  W  (lit)  100 


#  Digits /Sequence  Period.  K»q 
Symbol  Modulation 

f  Digit s/SyinboiaKBm 

#  Symbols/ Frame.  Kja 

f  Sequence  Periods/  Frame,  Kjm 


Periodic,  pseud  ■  random 


Variable  Symbol  Modutatim 


(C)  The  M5  experiment  had  two  parts,  each  hav¬ 
ing  distinctly  different  purposes.  The  objective  of 
the  M5A  experiment  was  to  measure  the  commu¬ 
nication  system  performance  at  low  sigual-to-noise 
ratios.  Consequently,  the  transmitter  power  was 
reduced  to  give  an  average  signul-to-noise  ratio  of 
0  db.  The  resulting  low  signal-to-noise  ratios 
increased  the  variance  on  the  transmission  meas¬ 
urements.  After  completion  of  the  M5A  portion, 
the  signal  power  was  increased  for  the  second  purl, 
of  the  M5  experiment  to  improve  the  quality  of 
the  transmission  measurements.  Because  of  the 
high  input  signal-to-noise  ratios  in  the  M5B 
experiment,  the  error  rate  was  reduced  to  a  point 
where  system  error  performance  could  not  be 
evaluated. 

(C)  The  M6  experiment  had  three  parts,  each 
representing  a  different  choice  of  system  param¬ 
eters.  In  the  M6A  experiment,  the  same  trans¬ 
mission  format  was  used  as  in  the  M5  experiment 
so  as  to  obtain  comparison  data  over  the  two 
distances.  The  sequence  duration  in  this  first 
experiment  was  known  to  be  too  short  for  reliable 
operation  at  the  42-nmi  range.  Further,  the 
average  signal-to-noise  ratios  at  Bimini  were  lower 
than  expected,  because  of  inadequate  source  drive 
at  Fowey  Rocks.  Consequently,  a  second  experi¬ 
ment,  M6B,  with  extremely  high  noise  and 
multipath  resistance  (and  low  data  rate)  was 
initiated.  Results  from  the  M6B  experiment 
indicated  that  its  design  was  too  conservative  and 
that  a  compromise  between  the  M6A  and  M6B 
parameters  was  in  order.  Subsequently,  the  third 
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M5A,  MSB,  M6A 

M6B 

1 - 

M6C 

Gross  Characteristics: 

Digit  Duration,  T  (ms) 

19 

67 

33 

Digit  Waveshape 

rectangular 

(wideband, 

coded) 

Sequence  Period,  (ms) 

286 

1000 

500 

Symbol  Duration,  Tgrn  (ms) 

57 

200 

100 

Frame  Component  Duration,  Tjm  (ms) 

1100 

4000 

2000 

Processor  Characteristics: 

A/D  Compression  Factor,  K 

Samples/Digit,  Nd_ 

Samples/Sequence  ’period,  NS(J 

Frames /Reference  Determination,  K 

8 

8 

4 

4 

14 

14 

60 

210 

210 

32 

16 

16 

Reference  Determination  Time,  (sec)  ‘ 

73 

128 

64 

SNR.n  Measurement  Averaging  Time  (min) 

8.  5 

14.  9 

7.  5 

Performance  Characteristics: 

Probe  processing  gain  (dB) 

30 

27 

27 

Symbol  processing  gain  (dB) 

17 

22 

19 

Theoretical  Minimum  SNHj  for  pe=.001  (dB) 
Transmission  Rale  (bit/sec) 

-7 

-12 

-9 

2.2 

.  63 

1.  25 

experiment,  M6C,  was  implemented  as  an  opti¬ 
mized  system  for  the  42-nmi  distance. 

(IT)  The  five  experiments  have  very  similar  trans¬ 
mission  formats;  indeed,  the  only  major  difference 
among  the  experiments  is  in  the  digit  structure 
and  duration.  Table  2  depicts  the  parameters  that 
were  common  to  all  experiments,  but  later  para¬ 
graphs  will  point  out  the  differences  among  the 
remaining  system  parameters.  Because  of  the 
common  features  of  the  experiments  and  the 
software  implementation  of  the  processor,  the 
transition  from  one  experiment  to  the  next,  was 
easily  accomplished. 

(U)  Although  five  different  experiments  were 
performed,  there  are  only  three  different  sets  of 
system  parameters,  with  the  M5A,  M5B,  and 
M6A  experiments  all  employing  the  same 
transmission  format.  Table  3  gives  the  significant 
differences  among  the  experiments. 

(C)  The  primary  difference  between  the  M5A, 
M5B,  and  MoA  format  and  the  M6B  and  M6C 
format  is  in  the  digit  durution  und  structure, 
with  a  digit  in  the  MOB  system  being  3.5  times 


longer  than  a  digit  in  the  M5A  system  and  a  digit 
in  the  M6C  system  being  1.8  times  longer  than  a 
digit  in  the  M5A  system.  The  purpose  of  the 
digit  lengthening  in  the  MOB  and  M6C  systems 
was  to  increase  the  symbol  processing  gains  and 
improve  the  performance  against  multipath. 

(C)  Of  the  last  two  systems,  MOB  and  MGC,  the 
MOB  has  the  highest  processing  gains  and,  cor¬ 
respondingly,  the  lowest  transmission  rate.  The 
MOO  is  exactly  twice  as  fast  and  has  8  db  less 
processing  gain  than  the  MOB  system.  Thus,  the 
MGO  system  can  be  considered  a  compromise 
between  the  short  digit  durations  of  the  M5A 
format  and  the  very  long  durations  of  the  MOB 
format. 

(U)  Both  the  MOB  and  the  M6C  systems  employed 
broadband  coding  of  the  digit  waveform  to  avoid 
the  reduction  of  symbol  bandwidth  with  the 
increase  in  symbol  duration.  This  was  necessary 
to  retain  the  advantages  of  a  broadband  trans¬ 
mission  against  noise  and  multipnth  when  the 
symbol  duration  is  long  relative  to  the  digit 
bandwidth.  A  comparison  of  the  digit  spectra  for 

CONFIDENTIAL 


.... ...  „  .  ..... 


.  s  ...  ... 


16 


WHMWWI' 


f 

J. 

I 


i 

t 


'  J'-trt-m*-  »i  ►, 


CONFIDENTIAL 


•T)  Digit  spectrum  for  M6C  tiptrtmtnt 


Figure  (>  ComimriMin  of  digit  spectra.  (Figure  unclassified.) 

ilu>  three  formats  is  shown  in  Figure  6.  Bemuse 
of  ?!i  •  tillering  of  the  source  and  reflector  and  the 
rec  .sei  analog  filters,  only  the  spectral  compon¬ 
ent.  in  a  100-hz  hand  about  the  carrier  are 
impoi  taut. 

U.  EXPEltIMEN'TAL  RESULTS 

tlT)  The  five  experiments  yielded  data  on  both 
acoustic  transmission  conditions  und  system  error 
performance.  In  this  set  lion,  only  measurements 
that  hear  on  the  evaluation  of  system  error 
performance  or  on  future  applications  on  the 
communication  system  are  discussed.  The  re¬ 
maining  acoustic  transmission  measurements  have 
been  made  available  to  other  scientists,  und  they 
may  be  reported  later. 

/.  Siijixil-lo-Xoisc  h ‘alio  Ui.il(»jrams 
<  U )  To  interpret  the  communication  performance 
data,  an  understanding  of  the  signal  and  noise 
environment  of  eaeh  experiment  is  necessary. 
Figures  7  and  8  depict  histograms  of  the  wideband 
received  -igmil-to-noise  ratio,  SNR,,,,  with  ordin¬ 
ates  of  both  frequency  in  percent  and  time  in 
hours.  The  ah-cissa  represents  SN R,„  averaged 
over  the  time  interval  shown  in  Table  d,  und  it 


Figure  7a.  Histogram  of  wideband  received  signul-to-noiso 
ratio:  Mo  A  experiment,  (Figure  unclassified.) 


Figure  7b.  Histogram  of  wideband  received  signal- to-noise 
ratio:  MoB  experiment.  (Figure  unclassified.) 


is  quantized  in  1-db-wide  bins.  Because  tho 
source  level  was  not  held  eoustunt  during  the 
experiments,  these  histograms  relied  only  the 
signal-to-noise  ratio  environment  as  seen  by  the 
receiver,  not  a  general  measurement  of  acoustic 
propagation.  Nevertheless,  they  indicate  the 
extent  of  data  available  ut  a  given  input  signal-to- 
noise  ratio,  so  the  statistical  significance  of  the 
error  measurements  can  be  understood. 

(C)  Figure  7  gives  the  signal-to-noise  ratio 
histograms  for  the  two  phase-  of  the  M5  series 
of  experiments.  As  mentioned  earlier,  the  source 
levels  for  the  M5A  experiment  were  deliberately 
reduced  to  achieve  the  low  average  signal-to- 
noise  ratios  shown  in  Fig.  7a.  Consequently,  a 
large  amount  of  data  with  reasonably  frequent 
errors  whs  obtained  in  the  M5A  experiment 
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Figure  811.  Hi-i«grum  of  wideband  received  signal-to-noise 
ratio:  MtiA  ex|x  riincnt .  (Figure  unclassified.) 


Figure  Hb.  ULstogram  of  wideband  received  xignal-to-noUe 
ratio:  MtiB  <-X|R'riiiH'iit.  (Figure  unclassified.) 


Figure  7I>  sixes  the  corresponding  histogram  for 
the  MoB  experiment,  in  which  the  source  wus 
driven  at  maxiitiuni  available  lex-els,  yielding 
about  a  12-db  increase  in  mean  SNR,,,  ox-er 
that  of  the  M5A  experiment.  Because  only 
1.2  percent  of  the  M5B  data  correspond  to 
SXR,„  beloxv  0  db,  the  error  performance  data 
for  the  M5B  experiment  xvere  not  statistically 
significant  and  xxere  not  processed  further.  Thus, 
the  MoB  experiment  does  not  appear  in  the 
subsequent  presentation  of  the  error  performance 
data. 

(U)  Figure  8  gives  the  SNR,„  histograms  for  the 
three  M6  series  experiments.  The  MtiA  experiment 
was  o|N*ruted  at  signal-to-noise  ratios  that  were, 
only  slightly  (6  db)  xvorse  than  those  of  the  M5A 
experiment,  so  an  adequate  measurement  of  sys¬ 
tem  performance  could  be  achieved.  Because  of 
failure  in  an  amplifier  driving  the  Fowey  Rocks 
source,  very  loxv  signal-to-noise  ratios  were  avail¬ 
able  for  the  M6B  experiment,  us  shown  in  Fig.  8b. 


Figure  8c.  Histogram  of  wideband  received  signal-to-noise 
ratio:  M6C  experiment.  (Figure  unclassified.) 

These  low  signal-to-noise  ratios  motivated  the 
choice  of  system  parameters  for  the  MOB  experi¬ 
ment.  Later,  a  new  amplifier  was  installed  at 
Fowey  Rocks,  xvhich  improved  the  signal-to-noise 
ratio  situation  at  Bimini,  and  the  MtiC  experiment 
was  initiated.  As  shown  in  Fig.  8c,  the  MfiC 
experiment  experienced  about  the  same  S.\R,„ 
environment  as  the  MtiA  experiment  had  experi¬ 
enced  earlier. 

2.  Communication  System  Performance 
(U)  Evaluation  of  communication  system  per¬ 
formance  is  a  difficult  problem  when  one  is 
confronted  with  the  varying  conditions  of  the 
acoustic  medium.  The  technique  described  here  is 
based  on  a  well-known  theoretical  model,  and  it 
offers  several  advantages  in  terms  of  clarity  and 
convenience. 

(U)  The  primary  data  for  the  performance  meas¬ 
urement  cnn-ist  of  input  signal-to-noise  ratio 
averages  and  the  cumulative  error  count  E— 
Ep+E,  taken  over  the  averaging  interval.  Table  4 
gives  the  averaging  times  and  corresponding 
number  of  bit  decisions  for  each  experiment. 
Because  the  number  of  bit  decisions  is  known,  u 
probability  of  error  for  each  data  pair  can  be 
computed.  A  natural  question  arises  as  to  the  best 
way  of  presenting  the  relationship  between  average 
input  signal-to-noise  ratio  and  probability  of  error 
when  thousands  of  such  pairs  arc  available. 

(U)  One  method  of  evaluating  the  system  per¬ 
formance  would  be  to  sum  all  of  the  error  counts, 
/?,  for  a  fixed  average  signal-to-noise  ratio  and 
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T\bif.  4.  Averaging  time*  of  SXIi and  bit  deei- 
gioim.  ( Table  cfa*»ified  Confidential.) 


MM,  MSI.  M6A 

M6I 

M6C 

SKR.  Av«r  Aging  t  lac  (Min) 

8.5 

14.9 

7.5 

No.  of  bit  doc 1* ions  in 
Averaging  tiao 

2205 

1085 

1085 

then  divide  bv  the  total  number  of  bit  derisions 
involved  to  obtain  a  probability  of  error.  The 
disadvantage  of  this  technique  is  that  it  does  not 
Indicate  the  variations  in  system  error  performance 
over  time  To  depict  such  variations,  the  error 
performance  of  each  block  is  evaluated,  with  the 
results  placed  in  histogram  form  to  indicate  the 
frequency  of  occurrence. 

(U)  A  histogram  of  error  probabilities  would 
satisfy  the  above  requirement;  however,  the  non- 
nnifomiity  of  the  probability  scale  makes  the 
histogram  difficult  to  understand.  To  overcome 
this  problem,  a  conversion  of  the  measured 
probability  of  error  to  a  more  uniform  axis  is 
helpful.  A  particularly  helpful  conversion  from  a 
systems  analysis  standpoint  is  to  go  from  proba¬ 
bility  of  error  to  output  signal-to-noise  ratio, 
which  can  be  presented  in  decibels. 

(U)  Several  alternatives  exist  for  converting  from 
orroi  probability  to  output  signal-to-noise  ratio, 
but  the  simplest  technique  is  to  use  a  known 
relation  from  an  existing  theoretical  model.  In  this 
paper,  the  conversion  of  probability  of  error  to 
output  signal-to-noise  ratio  is  based  on  the  classical 
results  for  the  likelihood  ratio  receiver  for  a  known 
linear,  time-invariant  channel  with  added  white 
Gaussian  noise.  For  this  receiver,  the  probability 
of  error,  />E,  is  given  by 

pE=4(-vS)  (5») 

where 

(e_,5  j)/\2tr  d*  (*°) 

and  D  is  the  output  signal-to-noise  ratio.  For  this 
situation,  the  conversion  of  error  probability  to 
output  signal-to-noise  ratio  can  be  accomplished 
by  solving  Eq.  (60)  for  Z>.  That  is 

D=\*-'(p*)V  <•»> 

where  4>"'(r)  is  the  inverse  of  the  function  ♦(«) 
defined  in  Eq.  (60).  In  the  idealized  situation  of 
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Figure  9u.  Ill- digram  of  values  of  O:  8NRi«— — 4  db. 
(Figure  classified  Confidential.) 


Figure  9b.  Histogram  of  values  of  D:  8NRi.—  6  db. 
(Figure  unclassified.) 


-4  -i  0  2  *  •  •  10  t  dB) 

Figure  9c.  Histogram  of  values  of  D:  8NRi.=  -8  db. 
(Figure  unclassified.) 
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Figure  fid.  Histogram  i>f  values  of  l>:  SNRi„=  — 10  db. 
( Figure  unclassified.) 


the  model,  D  is  in  fact  the  mcan-squared-to- 
variitnre  ratio  of  the  inutehed-filtcr  output.  The 
motivation  for  the  above  choice  of  theoretieal 
model  comes  front  both  its  agreement  with  the 
assumed  conditions  of  the  real  acoustic  channel  and 
its  ease  of  use. 

(U)  The  analysis  technique  from  this  point  is 
straightforward.  The  error  probability  for  each 
block  of  data  is  converted  to  output  signal-to-noise 
ratio  I)  in  decibels  through  Eq.  (61).  Then, 
histograms  of  the  resulting  values  of  U  are  plotted 
for  each  signal-to-noise  ratio.  This  provides  a 
smooth  atid  easily  understood  histogram,  as  shown 
in  Fig.  9. 

(U)  One  disadvantage  of  the  histogram  of  output 
signal-to-noise  ratios,  I),  is  that  the  smallest 


Figure  10.  Mean  D  versus  SNH,„  (Figure  classified 
Confidential.) 


Figure.  It.  Standard  deviation  of  0  versus  SNR,,,.  (Figure 
classified  Confidential.) 

measurable  probability  of  error  is  ubout  0.001, 
which  results  from  the  fact  that  only  one  or  two 
thousand  bit  decisions  are  made  in  each  block. 
Such  a  probability  of  error  corresponds  to  a  value 
of  D  of  about  1 1  db,  so  the  histograms  arc  neces¬ 
sarily  truncated  at  that  value. 

(U)  Figure  10  depicts  the  means  of  the  output 
signal-to-noisc  histograms  for  the  four  communi- 
eution  experiments  (M5B  is  omitted,  because  of 
its  excessive  signal-to-noise  ratios).  Points  are 
shown  on  this  plot  only  if  at  least  10,000  bit 
decisions  were  made  at  the  given  SNRm-  Conse¬ 
quently,  each  point  shown  is  a  svatistiudly  sig¬ 
nificant  measure  of  the  system  performance. 
Figure  1 1  depicts  the  standard  deviation*  of  the 
same  histograms. 

(C)  Note  that  there  is  a  distinctive  leveling  off  of 
output  signal-to-noise  ratio  with  increasing  input 
signal-to-noise  ratio.  This  is  most  likely  to 
the  non-Gaussian  nature  of  the  noise,  which 
becomes  significant  at  higher  signal-to-noisc  ratios. 
Also  worthy  of  note  is  the  fact  that  probabilities  of 
error  of  0.01  are  achieved  at  input  signal-to-noise 
ratios  well  below  0  db.  Thus,  a  satisfactory  system 
can  be  made  to  operate  with  SNRm  below  0  db. 
The  leveling  off  of  the  M6B  curve  is  probably  due 
to  the  very  long  reference  determination  time 
involved  (128  sec). 

(C)  Another  way  of  interpreting  communication 
svstem  performance  is  to  compare  the  experi¬ 
mental  system  performance  with  that  of  the 
idealized  system  described  previously.  For  the 
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idea!  system,  the  output  signul-to-noise  ratio  D  is 
given  by 

Dm*~VVTJCm'Km£S*m  (62) 

At  a  fixed  SNR,„,  the  difference,  between 
UMr,,  and  l>  can  be  found  to  indicate  relative  sys¬ 
tem  performance: 


This  has  the  advantage  of  placing  systems  with 
different  transmission  rates  on  the  same  basis. 
For  example,  the  M6B  system  is  better  than  the 
M6C  system  in  absolute  performance,  because  of 
its  slower  transmission  rate.  When  compared 
relative  to  the  ideal  system  with  the  same  rate, 
however,  the  relative  performance  of  the  M6C 
system  is  better.  This  is  shown  in  Fig.  12  for  the 
four  systems. 

(<_’)  Of  the  four  systems  shown  in  Fig.  1 1  the  M5A 
system  has  the  best  relative  performance,  while 
M6B  has  the  worst  relative  performance.  This  is 
probably  because  of  the  increased  stability  of  the 
7-nitii  path  with  a  73-sec  probe  measurement 
nterval  over  the  42-nmi  path  with  a  128-soc 
probe  measurement  interval.  Similarly,  the  relative 
performance  of  the  MCA  and  M6C'  systems  is 
worse  than  that  of  the  N15A  system  strictly 
because  of  range.  The  relative  performance  dif¬ 
ference  between  the  M5A  and  M6C  systems 
shown  in  Fig.  12  is  not  considered  to  be 
significant. 

(C)  In  summary,  the  performance  of  the  communi¬ 
cation  systems  was  very  well  established  with 
large  amounts  of  data.  With  careful  choices  of 
system  parameters,  a  coherent  matched-filter 
communication  system  such  as  M6C  can  yield 
error  probabilities  of  0.01  at  input  signul-to-noise 
ratios  of  —10  db.  Such  a  system  would  be  operat¬ 
ing  with  a  performance  only  6  db  worse  than  that 
of  an  idealized  system. 

3.  Channel  Stability  Measurement * 

(C)  The  utility  of  the  matched-filter  technique 
studied  here  is  highly  dependent  on  the  stability 
of  the  medium.  If  the  acoustic  channel  changes 
significantly  from  the  time  of  the  probe  measure¬ 
ment  to  the  time  of  the  information  component 
processing,  then  the  filter  will  not  be  properly 
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Figure  12.  Mean  *  versus  SNR...  (Figure  classified 
Confidential.) 

matched  and  an  increase  in  errors  can  be  expected. 
On  the  other  hand,  if  the  acoustic  medium  is  very 
stable,  then  the  integration  time  of  the  probe 
component  can  be  increased  to  improve  the 
measurement  of  the  channel  digit  response. 

(U)  Because  of  the  importance  of  the  channel 
stability  to  future  applications,  the  five  experi¬ 
ments  included  an  elementary  measurement  of 
the  wideband  channel  stability,  specifically,  the 
normalized  cross  correlation  coefficient  p(k)  tie- 
scribed  in  Section  111b.  Although  correlation 
coefficient  data  are  available  for  the  full  duration 
of  each  experiment,  only  a  few  days  of  such  data 
for  each  experiment  have  been  analyzed.  These 
days  were  chosen  from  continuous,  high-quality 
data  to  determine  any  daily  fluctuations.  Only 
data  with  SNK*,,  greater  thun  10  db  were 
considered. 

(U)  The  most  convenient  way  to  understand  the 
variations  of  the  acoustic  channel  is  to  construct 
histograms  of  the  correlation  coefficient.  Figure  13 
depicts  histograms  of  p(k)  for  the  NI5A,  M6A,  and 
M6B  experiments.  The  histograms  represent 
one-day  intervals  (about  1,440  correlation  values 
per  histogram)  with  a  bin  size  of  0.01.  Note  that 
the  abscissa  on  each  plot  begins  at  0.5,  not  zero. 
Since  M5A  and  M6A  used  the  same  transmission 
parameters,  the  decrease  in  mean  value  can  be 
attributed  to  the  increased  distance,  42  nmi 
instead  of  7  nmi. 

(C)  Figures  13b  and  13c  depict  the  p(k)  histograms 
for  the  M6A  and  M6B  experiments.  The  obvious 
widening  of  the  M6B  histograms  is  due  to  the 
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c)  M6B  Experiment 


Figure  13.  Histogram*  of  pit).  (Figure  unclassified.) 

increased  measurement  interval  durations,  128 
sec  instead  of  73  sec.  The  lower  correlation  values 
of  the  MOB  experiment  indicate  that  the  measure¬ 
ment  of  the  digit  response  is  less  accurate  in  the 
MOB  system  than  in  the  MOA  system.  This  offers 
at  least  a  partial  explanation  for  the  poor  relative 
performance  of  the  MOB  system  shown  in  figure 
12.  The  correlation  function  histograms  for  the 
M5A  and  M5B  experiments  are  nearly  identical, 
as  the  distance  and  integration  times  arc  similar 
for  the  two  experiments.  Also,  the  histograms  for 
the  M6A  and  M6C  experiments  are  similar. 

V.  SPATIAL  STABILITY  MEASUREMENTS 

(U)  The  preceding  sections  described  the  opera¬ 
tion  ami  |K>rformunce  of  an  underwater  acoustic 
communication  system  operating  between  two 


fixed  points.  In  a  practical  submarine  communica¬ 
tion  system,  one  or  both  ends  of  the  acoustic 
channel  is  in  motion,  so  spatial  variations  as  well 
as  temporal  variations  are  important.  This  sections 
describes  tw  o  measurements  of  the  spatial  stability 
of  the  medium  and  provides  the  foundation  for 
the  extension  of  the  fixed-site  communication 
technique-.:  previously  discussed  to  the  submarine 
communication  problem. 

(U)  The  fixed-site  experiments  have  shown  the 
temporal  stability  of  the  medium  to  be  sufficient 
to  allow  integration  limes  of  the  order  of  one 
minute.  From  a  purely  geometric  point  of  view, 
the  medium  should  also  be  stable  under  spatial 
displacements  that  are  common  to  submarine 
platforms.  For  example,  a  submarine  on  a  10-knot, 
zero  Doppler  track  at  a  range  of  100  nnti  subtends 
less  than  0.1  deg  of  arc  in  one  minute.  The  wide¬ 
band  characteristics  of  the  acoustic  channel 
would  not  be  expected  to  change  significantly 
under  such  a  displacement.  Nevertheless,  careful 
measurements  have  been  mude  to  validate  this 
intuitive  understanding. 

(U)  Two  measurement  techniques  were  employed 
to  study  the  spatial  stability  of  the  acoustic 
channel.  The  first  used  coherent  cross  correlation 
of  received  waveforms  without  any  coherent 
integration,  and  it  yielded  anomalous  results.  The 
second  technique  used  coherent  cross  correlation 
with  13.4  db  of  coherent  integration,  and  this 
produced  consistent  results.  In  the  following  sec¬ 
tions,  each  measurement  technique  is  described 
in  detail,  anti  an  explanation  of  the  differences  in 
their  results  is  offered. 

A.  COHERENT  moss  C  ORRELATION  WITHOUT 
COHERENT  INTEGRATION 

(U)  In  this  measurement  technique,  the  spatial 
stability  of  the  medium  was  measured  in  terms  of 
normalized  correlation  coefficients,  pi*(i),  analo¬ 
gous  to  the  correlation  coefficient,  p(.k),  studied  in 
the  M5  and  M6  experiments.  Because  of  the 
presence  of  the  Doppler  effect,  however,  the 
definition  had  to  be  modified  slightly. 

(U)  The  transmission  was  a  15-digit,  pseudo¬ 
random  sequence  similar  to  that  of  the  M5A, 
M5B,  and  MtiA  experiments.  The  sequence  was 
scaled  from  the  previous  carrier  frequency  of  420 
hz  to  a  center  frequency  of  350  hz,  to  be  com- 
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putible  witli  tin*  acoustic  source.  The  resulting 
sequence  has  a  period  Tm  of  M  un<1  “ 

zero-to-zero  bandwidth  of  87.5  hz.  and  it  was 
trail,  milted  eontinuously  from  the  lowed  souree 
during  each  experiment. 

(U)  Is't  /,(/)  be  the  received  signal,  bandpass- 
Kltered  in  a  100-liz  band  centered  on  the  carrier 
frequence  of  ::50  hz.  Consider  a  segment  of  r,(<) 
of  duration  7',,  that  contains  an  integer  number  of 
sequence  periods  at  zero  Doppler,  that  is,  T„— 
l-T ^  for  some  integer  k.  The  normalized  cross 
correlation  function  p,tr)  of  this  segment  with  the 
input  r,(t)  is  given  by 


small,  the  consequences  of  such  time-axis  distor¬ 
tions  are  significant  when  coherent  techniques  are 
employed. 

(U)  If  the  velocity  of  the  source  and  the  average 
sound  velocity  were  known,  the  normalized 
correlation  function  of  Eq.  (64)  could  be  evaluated 
at  multiples  of  the  Doppler-distorted  sequence 
duration,  TV  Unfortunately,  neither  variable  is 
sufficientlv  well  determined  to  calculate  T'..  ade¬ 
quately.  In  the  September  1978  experiments,  a  new 
correlation  coefficient,  p,*(0,  was  defined  so  as  to 
alleviate  the  effects  of  Doppler.  Let, 

Pl*(0=max|p1|iT,(l-<i)ll;  -/<«<+*  (66) 


rt<.+rH 

Jit 


i,(tU,*(t+T)dt 


h+T„ 


|/i(OIL’ 


/*#o+T 

!(lt 
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l#+7'»a  +  r 


\MtWdt 


(64) 


where  tn  is  the  starting  time  of  the  selected  seg¬ 
ment.  Note  that,  in  general,  p,(t)  is  complex 
valued.  H  there  is  no  Doppler  and  the  medium  is 
perfectly  stable,  then  p,(r)  will  be  periodic,  with 
the  same  period  as  the  sequence,  and  it  will  attain 
a  maximum  value  of  unity.  If  there  is  no  Doppler, 
but  the  medium  is  changing,  values  of  p,(r)  taken 
at  multiples  of  the  sequence  period  are  a  measure 
of  the  stability  of  the  channel.  Consequently,  the 
correlation  function  evaluated  at  multiples  of  the 
segment  duration  p,(iTv),  «=1,  2  .  .  .  ,  will  also 
be  an  indication  of  channel  stability. 

(U)  Since  /,(()  was  generated  by  a  moving  source, 
the  effects  of  Doppler  will  be  present,  even  on  a 
nominal  “zero  Doppler”  track.  A  convenient 
approach  to  understanding  the  effects  of  Doppler 
is  to  consider  it  as  producing  either  a  time  com¬ 
pression  or  expansion  of  the  received  signal.  Let 
7'„  be  the  duration  of  a  segment  of  x,(t)  at  zero 
Doppler  and  let  T'„  be  the  duration  of  the  same 
segment  as  received  from  the  source  approaching 
the  hydrophone  ut  a  velocity  r.  Ihen 

T'm=(\-v/c)T„  (65) 

where  c  is  the  average  speed  of  sound  in  the  same 
units  us  r.  If  the  source  is  approaching  the  receiver 
at  10  knots,  7*  is  reduced  by  0.84  percent. 
Similarly,  if  the  source  is  leaving  the  receiver  at 
10  knots,  Tm  will  be  increased  by  0.84  percent. 
Although  the  percent  variation  may  appear  to  be 


where  i=l  ...  7. 


(67) 


v*  is  the  largest  anticipated  Doppler  rate  and 
c,„  is  an  estimate  of  the  average  speed  of 
sound.  The  purpose  of  selecting  the  maximum 
magnitude  of  pt(r)  in  an  interval  is  to  ensure  that 
the  correlation  coefficient  corresponds  to  the 
correlation  value  at  correct  time  alignment  in 
the  presence  of  Doppler. 

(U)  Two  important  differences  between  the  cor¬ 
relation  coefficient  p(k)  of  the  M5  and  M6  experi¬ 
ments  and  the  correlation  coefficients  Pi  *(.»)  of  the 
September  1974  experiment  must  be  noted.  First, 
in  the  M5  and  \16  experiments,  p(k)  was  derived 
from  a  cross  correlation  of  two  coherently  inte¬ 
grated  waveforms  p(k,i),  p(k—\,>)  in  Eq.  (36) 
[whereas  p,*(»)  is  derived  from  a  cross  correlation 
of  two  unintegruted  input  wuveforins,  /,(<)> 
in  Eq.  (64)[.  Consequently,  the  coefficients 
p,*(i.)  are  much  more  sensitive  to  both  noise  and 
short-term  channel  variations  than  the  coefficient 
p(k)  measured  previously. 

(U)  The  second  important  difference  between  the 
measurement  of  p,*(<)  in  the  September  1974 
experiment  and  the  measurement  of  p(k)  in  the 
MS  and  M6  experiments  is  that  pi  *(»)  is  insensitive 
to  pure  time  shifts  in  the  medium.  That  is,  p,  (i) 
can  be  close  to  unity  for  a  medium  exhibiting 
rapid  changes  in  time  delay  because  of  the  maxi¬ 
mization  included  in  the  definition  of  p,  (i). 
Fortunately,  many  previous  transmission  measure¬ 
ments  have  shown  that  such  changes  do  not  exist 
in  the  ocean  acoustic  medium. 


CONFIDENTIAL 


CONFIDENTIAL 


l-OWKY  ROCKS 

O 


R  liK |;R  KXC  K:  C  AC.S  1248 
DEPTHS  IN  FEET 


Figure*  14.  Vessel  track  for  the  September  1973  experiments.  (Figure  unclassified.) 
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Figure  ISo.  Average  correlation  coefficients  versus  spatial 
displacement,  Sept.  14,  1073.  (Figure  unclassified.) 


Figure  tub.  Average  correlation  coefficients  versus  spatial 
displacement,  Sept.  21,  1973.  (Figure  unclassified.) 

(U)  The  preliminary  spatial  stability  measure¬ 
ments  were  conducted  with  a  towed  HX90  source 
in  the  Straits  of  Florida  during  September  1973. 
Figure  14  depicts  the  vessel  track  during  these 
experiments.  The  source  was  towed  at  approxi¬ 
mately  200-ft  depth  with  vessel  speeds  of  2.5,  5, 
und  10  knots.  Signals  from  the  source  were 
received  at  a  hydrophone  approximately  7  nmi 
from  the  Fowey  Rocks  Light  and  were  sent  via 
cables  to  the  laboratory  on  shore. 

(U)  Figure  15  depicts  the  average  correlation 
coefficients  obtained  during  three  experiments  in 
September  1973.  Here  Tm=5  sec  and  i  ranges 
from  1  to  7,  corresponding  to  a  40-sec  maximum 
extent.  The  coefficients  p,*(i)  are  plotted  as  a 
function  of  distance,  not  time,  so  as  to  indicate 
the  stability  for  a  fixed  spatial  displacement, 

CONFIDENTIAL 


Figure  13c.  Average  correlation  coefficients  versus  spatial 
displacement,  Sept.  23,  1973.  (Figure  unclassified.) 
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Figure  19.  Average*  of  average  correlation  coefficient* 
versus  spatial  displacement.  (Figure  unclassified.) 

independent  of  source  velocity.  Because  only  40 
sec  of  data  could  be  processed,  the  spatial  varia¬ 
tions  are  measured  over  shorter  displacements  for 
lower  source  velocities.  Data  for  which  Pi*(i)<[0.8 
have  been  eliminated  from  the  plot,  as  such  a  value 
is  indicative  of  an  inadequate  signal-to-noise  ratio. 

(U)  The  data  given  in  Fig.  15  exhibit  several 
obvious  anomalies.  First,  the  values  of  the  cor¬ 
relation  coefficients  do  not  agree  for  the  same 
spatial  displacement.  If  these  correlation  coeffi¬ 
cients  actually  measured  the  spatial  variations  of 
the  medium,  coefficients  corresponding  to  the 
same  spatial  displacement  should  be  equal,  which 
is  not  the  case.  Second,  on  different  days  the 
relative  positions  of  the  curves  for  different  source 
speeds  change.  For  example,  on  September  14,  the 
2.5-knot  data  were  the  most  stable,  whereas  on 
September  21  and  25  the  2.5-knot  data  were  the 
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Figun  17.  Averages  of  average  correlation  coefficients 
versus  time  displacement.  (Figure  unclassified.) 

•east  stable.  Similarly,  the  relative  positions  of  the 
5-  ami  10-knot  rims  change  from  one  day  to  the 
other. 

(U)  Figure*  1G  gives  the  average  of  the  results 
depicted  in  Figs.  15a.  15b,  and  15c.  As  in  the  one- 
day  plots,  the  correlation  coefficients  do  not  appear 
to  depend  on  spatial  displacement.  Further,  the 
correlation  coefficients  are  smaller  for  low  source 
speeds  than  for  higher  source  speeds.  Neither 
result  agrees  with  common  sense. 

(U)  In  Fig.  17  the  average  correlation  coefficients 
of  Fig.  Iti  are  plotted  as  a  function  of  time  displace¬ 
ment  instead  of  spatial  displacement.  The  close 
agreement  lie! ween  the  curves  for  different  source 
speeds  is  apparent.  This  indicates  that  the  cor¬ 
relation  coefficients  p, *(/)  are  dependent  on  time 
displacement ,  not  spatial  displacement  as  orig¬ 
inally  intended.  Also,  the  time  stability  is  dras¬ 
tically  reduced  from  the  values  obtained  in  the 
extensive  fixed-site  experiments  (M5  and  Mb), 
when*  correlation  coefficients  of  0.9  were  obtained 
over  longer  intervals  (7H  and  64  sec.). 

(U)  Because  of  the  above  anomalies  in  the  data, 
a  search  for  possible  explanations  was  initiated.  A 
brief  study  was  made  of  a  cw  transmission  from 
the  towed  source  to  determine  any  peculiarities  in 
the  source  motion.  I'll  is  study  showed  that  varia¬ 
tions  of  the  source  about  its  intended  linear  track 
were  insignificant. 

(U)  Flu*  most  probable  causes  of  the  anomalies 
in  the  data  are  noise  and  surface  modulation. 
Although  the  signal-lo-uoi.se  ratio  during  eacli 


experiment  was  estimated  to  la-  in  excess  of  10  db, 
no  quantitative  measurement  such  as  those  in  the 
M3  and  MG  experiments  was  available.  Because 
no  coherent  integration  was  used,  the  correlation 
coefficients  p,*(<)  are  particularly  sensitive  to 
noise.  Although  previous  measurements  in  the 
Straits  have  indicated  that  the  extent  of  surface 
modulation  is  insufficient  by  itself  to  account  for 
the  observed  results,  surface  modulation  would 
hIso  reduce  the  values  of  the  correlation  coeffi¬ 
cients.  Both  of  these  effects  are  independent  of 
spatial  displacement,  and  they  vary  in  extent  on  a 
dny-by-day  basis,  which  is  in  agreement  with  the 
data  obtained. 

(U)  Although  the  experiments  described  in  this 
section  were  unsuccessful  in  measuring  the  spatial 
stability  of  the  medium,  the  lesson  learned  is 
important.  That  lesson  shows  the  importance  of 
substantial  coherent  integration  in  ihe  measure¬ 
ment  of  medium  stability,  so  the  effects  of  nniso 
and  surface  modulation  do  not  influence  t lie 
stability  measurement.  As  a  result  of  the  Septem¬ 
ber  1971  experiment,  a  second  measurement, 
technique  employing  coherent  integration  was 
developed,  and  this  yielded  a  valid  measurement 
of  spatial  stability. 

II.  COHERENT  CROSS  CORRELATION  WITH  COHERENT 
INTEGRATION 

(U)  After  the  preliminary  spatial  stability  meas¬ 
urements  in  September  197K,  development  of  a 
measurement  technique  including  coherent  inte¬ 
gration  was  initiated.  The  requirement  for  coherent 
integration  was  based  on  two  considerations. 
First,  the  September  1 973  experiments  show  ed  that 
noise  and  surface  modulation  precluded  measure¬ 
ment  of  the  spatial  stability  when  no  coherent 
integration  was  used;  and  second,  an  opportunity 
to  study  spatial  stability  at  long  ranges  in  the  deep 
ocean  was  available.  The  input  signal-to-noise 
ratio  at  these  ranges  required  c  oherent  integration 
to  obtain  adequate  representation  of  the  signal. 
The  resulting  measurement  technique,  which 
employs  coherent  cross  correlation  with  coherent 
integration,  is  described  below  . 

(U)  Because  of  the  need  to  perform  coherent 
integration  as  well  as  to  accommodate  Doppler 
effects,  fast  Fourier  transform  techniques  (eft) 
were  employed.  These  techniques  provided  two 
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constraints  on  the*  measurement  process.  First,  the 
coherent  integration  time  was  limited  to  about  -50 
see  because  of  the  memory  limitations  of  the 
processing  equipment.  Secotul.  the  transmission 
format  was  selected  to  have  an  integer  number  of 
sequence  periods  contained  in  the  transform 
interval,  so  that  a  close  alignment  of  the  signal 
and  fkt  s|M'(  tral  lines  would  be  obtained.  Once 
these  two  constraints  were  satisfied,  the  applica¬ 
tion  of  fkt  techniques  allowed  straightforward 
measurement  of  both  the  spatiul  stability  and  the 
input  signal-to-uoise  ratio. 

(U)  The  transmission  for  the  e\|H>riment  was  a 
3 1 -digit  linear,  maximal  sequenee  with  15  carrier 
cycles  per  digit.  Modulation  of  the  350-hz  carrier 
was  accomplished  by  complementary  phase  tech¬ 
niques.  so  as  to  yield  a  high-level  carrier  and  yet 
retain  a  higher  total  power  than  with  amplitude 
modulation.  The  transmission  is  given  by 

m(l )—  Xj  i  «(0  cos  2»/of+ll  - fl(01  sm  2r/0< } 

X  {U{iT)-U[(i+\)T]\  (68) 

where  /„  =  dot)  lix  is  the  carrier  frequency, 
7’=h.042s57  see  is  the  digit  duration.  |a(/)l«(0« 
|0,ll  |  are  the  am  sequence  coefficients,  and  U{t) 
is  the  unit  step  function.  This  transmission  has 
the  same  power  spectrum  as  an  amplitude- 
modulated  sequence  of  the  same  structure,  but 
it  has  approximately  :t  db  higher  total  power. 
Because  /„(/)  is  periodic,  with  a  period  of  1. 329567 
see.  its  spectrum  has  a  line  spacing  of  0.752125  lix. 
The  transmission  has  a  zero-to-zero  bandwidth 
of  46.667  hz,  determined  by  the  digit  duration,  T. 

(II)  After  initial  analog  bundpass  filtering  to  a 
50-hz  bund  about  350  hz,  the  received  signal 
y, (/)  is  sampled  at  a  rate  of  4 fnc  where  ftrt  is  the 
receiver  center  frequency  (  frrr--i 50  hz  for  zeto 
Doppler).  The  resulting  samples  were  compressed 
in  the  manner  defined  in  Eqs.  (26)  and  (27), 
but  with  /CP=10.  This  reduces  the  number  of 
samples  to  agree  with  the  sampling  theorem 
requirements  and  yields  186  compressed  sumples 
per  sequence  period.  For  convenience  in  the 
subsequent  dismission,  these  samples  are  consid¬ 
ered  in  terms  of  2.048  long,  complex  vectors, 
>Z(k,  «)!<-  0  .  .  •  2,047),  where  the  k  indicates 
the  fctli  transform  interval.  That  is,  the  sample 
following  Z(k,  2,047),  is  Z(k+ 1,  0).  and  so  forth. 


(U)  An  input  vector  of  2,048  complex  samples 
Z(k,  <)  contains  approximately  22  sequence  peri¬ 
ods  of  the  transmission  at  zero  Doppler.  Specifi¬ 
cally,  the  2,048  complex  samples  correspond  to  a 
29.257143-sec  time  interval,  whereas  22  sequence 
periods  occupy  29.228474  sec,  a  0.1  percent 
difference.  Consequently,  the  spectrul  lines  of  a 
2,048-complex-point  fft  will  very  nearly  align 
with  the  signal  spectrul  lines.  The  actual  offset 
between  the  7th  signal  line  from  the  carrier  and 
the  cot  responding  transform  line,  &Trr(')  is  given 
by 

(1  22  \ 

1.329567  ~ 29.257143/ 1  (69) 

=0.000728/  hertz  (70) 

For  the  energetic  lines  near  the  carrier,  this  offset 
is  less  than  half  the  fft  line  spacing  of  0.034180  hz, 
and  it  will  be  neglected  in  the  subsequent  dis¬ 
cussion. 

/.  Ope  rat  inn  at  Zero  Doppler 

(U)  To  understand  the  operation  of  the  measure¬ 
ment  technique,  first  consider  its  operation  with 
zero  Doppler.  Because  the  transform  interval 
contains  nearly  22  sequence  periods,  every  22nd 
transform  spectral  line  from  the  carrier  will 
contain  signal  energy.  The  intervening  21  trans¬ 
form  lines  will  contain  only  noise  energy.  Conse¬ 
quently,  a  processing  gain  of  22  (13.4  db)  can  be 
achieved  by  considering  only  every  22nd  trans¬ 
form  spectral  line  about  the  carrier. 

(U)  Let  |Z(*. /)|- 1,024 <(<1,0231  be  the  dis¬ 
crete  Fourier  transform  of  the  input  vector 
|z(J(r,  r)|i=0  .  .  .  2,0471  with  the  i  index  on 
\Z(k,  i)  |  selected  so  the  Z(k,  0)  corresponds  to  the 
carrier  line  at  zero  Doppler.  Define  a  new  spectrum 
{S(*,  ()|— 32<t<31 1  from  \Z(k,  i)j  by  selecting 
every  22nd  spectral  line  of  \Z(k,  01  follows: 

S(k,  i)=Z(k,22  i);  -32<»<31  (71) 

From  the  discussion  of  the  preceding  paragraph, 
\S(k,  /)!,  will  contain  signal  energy  with  a  13.4-db 
improvement  in  signal-to-noise  ratio  over  that  of 
\Z(k,  01,  because  of  the  elimination  of  21  lines 
containing  noise  alone. 

(U)  The  measurement  of  spatial  stability  is 
accomplished  by  calculating  the  normalized  cross 
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correlation  viwf"  'it.  p.*(i),  between  the  signal 
waveform  in  the  X'l li  transform  inlervat  with  that 
of  the  ik~  Dili  transform  interval.  I^ct  «((*,  /), 
#==0  .  .  .  127,  he  the  inverse  iliserete  Fourier 
transform  of  the  X't l»  -■iiriinl  speetrnm  S(k,  i).  The 
normali/.eil  < loss  correlation  function  p(k,j)  be¬ 
tween  *{k.  i )  and  x(k  I,  <)  is  defined  by 

*  Vvt  A* —  1 ,  i  +  j) 

- ,,  —  (72) 

1.  »)i* 

V  II 

where  the  (<+./)  index  is  taken  nimlnlo  128  so  that 
pik.j)  is  the  result  of  a  circular  correlation  Be¬ 
cause  the  transform  length  is  not  an  exact  multiple 
of  the  .sequence  period,  p(k, 0)  does  not  give  the 
correlation  corresponding  to  exact  alignment 
between  tin*  two  waveforms.  Consequently,  a 
maximization  process  must  be  performed  "to  deter¬ 
mine  the  correlation  coefficient,  ?i*(k),  which 
represents  the  medium  spatial  stability: 

pi*(k)  —  max  p(k.j)  ;  «<j<63  (7:?) 

The  coefficient  *►_,*(/-)  has  the  same  ambiguity  in 
terms  of  pure  lime  shifts  that  thcpi*<>)  coeffic  ients 
had  in  the  preliminary  experiments.  As  before, 
this  ambiguity  is  known  to  be  insignificant,  on  the 
basis  of  other  experimental  data. 

(U)  To  properly  interpret  the  correlation  coeffi¬ 
cient,  ?i*(k),  a  quantitative  knowledge  of  the 
signal-to-noise  ratio  of  both  the  S(k,  i)  and 
S(k—  I.  i )  vectors  is  needed.  The  signal  power 
SPOt)  can  be  computed  directly  from  the  signal 
spectrum  S(k,i). 

SP (*)-  £)  StA-,  »V  (74) 

<  ™  —  a  j 

To  determine  the  noise  power,  the  power  in  n 
transform  spectral  line  in  between  the  signal  lines 
is  calculated. 

XP(X)  £  \Z(k,22i-^U)\i  (75) 

i »  —  i j 

If  it  i>  assumed  that  the  noise  is  white  (at  least 
relative  to  the  signed  line  spacing  of  0.752691  hzl, 
SV(k)  gives  a  measure  of  the  noise  content  of 
S(k.  <).  Consequently,  the  signal-to-noise  ratio 
S.\R, lUtc/-)  of  .S’ (Ax  /)  is  simply  the  ratio  of  SP(j t) 
to  NP(X-) 


SXRlltt,<  k) 


SPcX) 

XP(X-) 


The  input  signal-to-noise  ratio  SX  R,„(Jt)  is  22 
times  (1.1.4  dh)  worse*  than  SN  H, „„(/(•). 

SXRm(t)  ■-=  ~  SN  H„„, (k)  (77) 

Thus,  l)oth  the  input  and  output  signal-to-nohc 
ratios  can  be  quantitatively  measured  in  each 
transform  interval 

2.  Operation  with  Xomero  Doppler 
(U)  Measurement  of  the  correlation  coefficient 
Pi*(k)  and  the  input  and  output  d^mil-to-noise 
ratios  is  not  significantly  more  difficult  when  the 
Doppler  is  nonzero.  Let/,/  be  the  Doppler-shifted 
carrier  frequency 

/  t  U  ,  re  \ 


If  the  difference  between  //  anti  the  receiver  center 
frequency,/,.,.,.,  is  small,  then  the  effect  of  Doppler 
on  the  received  signal  is  approximately  a  frequency 
translation.  Such  a  frequency  translation  nets  as 
u  fixed  offset  between  the  zero- Doppler  location 
of  the  signal  spectral  lines  and  the  actual  retched 
spectral  lines. 

(U)  Consider  the  situation  where/,/  and/rrr  differ 
by  At  It/.. 

A/=/c/— /r«  (76; 

Then  if  the  effect-  of  Doppler  are  considered  ns  a 
frequency  translation,  the  signal  spectrum  will  be 
displaced  by  /•'  tniii-form  spectral  lines. 


T>  6.034180  ^80) 

In  general,  the  displacement  r'  will  not  be  an 
integer.  The  integer  offset.  <*,  to  the  cloxmt 
transform  spectral  lino  can  be  found  from 

r  *  =  [<•' +  0.5)  (81) 

where  the  square  brackets  |/|  indicate  the  greatest 
integer  part  of  /.  In  subsequent  discussion,  r*  is 
milled  the  vernier  frequency  index. 

(lT)  (liven  the  vernier  frequency  index,  r*,  Eq. 
(71)  for  the  signal  spectrum  can  be  modified  to 
incorporate  the  effects  of  Doppler. 
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Fi|>nri'  IS.  Vessel  track  fur  the 
January  1 1*74  experiment.  (Figure 
uuclu—ificd.) 


Tru«  North 


T hi  ci  rainon*  of  fhi»  trick  cannot  b«  •Iran 
iu«  to  aacwrlty  ciMaiiaratiam 


Slk,  i)  =  ZAk.  22/+/'*)  (82) 

Suuiliirl \ .  E(|.  (75)  for  llie  noise  power  heroines 

Nl'U)  S  22/4-11+/*)^  (83) 

l=-3.' 

Willi  ilie-e  two  modification'*  of  the  zero- Doppler 
eipml  ioiis.  the  derivation  of  I  lit*  correlation 
roeHiricht  p_. *( A')  and  the  Mgnul-lo-nuise  ratios 
remain  \alid  in  tin*  presence  of  Doppler. 

( l’)  The  -ignilirnnt  remainin':  problem  is  the  deter- 
inimilion  of  the  vernier  frequency  index,  r*,  by 
the  measurement  sv-tem.  One  approach  is  to 
search  the  received  signal  spectrum  Z(h\  i)  for  the 
largest  -ignal  line  (presumably  the  carrier  at 
l»  i  s.)  and  calculate  If  and  /•*.  Unfortunately, 
'elect i\ e  fjitlinir  due  to  multipath  will  occasionally 
reduce  the  carrier  level  below  that  of  the  adjacent 
sidebands  and  the  system  will  obtain  an  incorrect 
value  for  i  *. 

lU)  To  avoid  this  problem,  the  measurement 
system  determined  r*  on  the  basis  of  broadband 
'iutial-to-noisc  ratio.  Specifically }  the  following 
quantity  i«  calculated 

ti3 

S  M  22/+r)|a 

SNK„U1  (lr,  /•)—  ,a!’u  -  (84) 

jZ/Ak,  2214^+11)1* 

<=u 

for  —  !()</  <  II.  The  expression  SNR„u,U\  r) 
given  in  K<|.  (84)  is  simply  the  output  signal-to- 
noi'C  ratio  if  r  were  tin1  correct  vernier  frequency 
index,  ri.c  range  of  /•  is  constrained  to  be  within 
one  signal-line  separation  (0.752091  Itz,  or  approx¬ 
imately  a  0.4-knot  Doppler  uncertainty),  as 
SN  K„,„ (/■ ■,  / )  will  be  periodic  because  of  the  regular 
spacing  of  the  signal  spectral  lines.  In  actual 
measurements,  the  receiver  frequency,  is 
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udjusted  to  be  within  about  0.1  hz  of  so  this 
constraint  in  r  is  not  a  problem. 

(U)  To  determine  the  value  of  r*,  the  processor 
selected  the  value  of  r  that  yielded  the  highest 
value  of  SNR„„,(fc,  r).  Because  SNR, m(k,  r )  is 
based  on  the  broadband  signal  spectrum  and  not 
on  any  single  spectral  line,  the  resulting  frequency 
offset,  /■*,  is  not  disturbed  by  selective  fuding. 
Given  /■*,  the  calculation  of  and  the  signnl- 

to-noise  ratios  is  easily  accomplished.  Indeed,  the 
measurement  system  performed  all  necessary 
calculations  in  real  time,  producing  new  values  of 
Pa*(it)  and  SNRoutU*)  every  29.257143  sec. 

3.  Experimental  llem/lte 

(U)  During  January  1974,  the  HX90  acoustic 
source  was  towed  in  the  deep  ocean  between 
Eleuthera  and  Bermuda.  The  approximate  vessel 
track  relative  to  the  hydrophone  is  shown  in 
Fig.  18.  Signals  from  the  source  were  received  at  a 
fixed  hydrophone  and  processed  with  the  coherent 
integra*io;,  technique  described  above.  The  ranges 
of  the  experiment  varied  from  0  tff400  nmi,  with 
most  of  the  experiment  being  conducted  at  ranges 
between  300  and  400  nmi.  During  the  experiment, 
the  vessel  speed  was  maintained  at  approximately 
6  knots. 

(U)  The  digital  processor  was  in  nearly  continuous 
operation  during  the  12  days  of  the  experiment, 
yielding  correlation  coefficients  #>j*(A*)  and  signal- 
to-noise  ratios  every  29.257143  sec.  These  and 
other  results  were  recorded  on  digital  magnetic 
tape  and  simultaneously  plotted  on  a  multipoint 
recorder.  Figure  19  depi.  ts  an  annotated  typical 
multipoint  plot  from  the  processor.  As  in  the  M5 
and  M6  experiments,  this  plot  was  an  important 
aid  in  maintaining  quality  control  during  the 
experiment. 
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Figure  22a  lli-lograni  nf  p;*ik)  values  for  full  Diipplit 
('■xiditioiis.  1  Figure  unclassified.) 


histograms  of  p.» * (/*) ,  much  us  in  the  M5  «n<l  M6 
experiment-.  Because  noise  cun  reduce  the  values 
of  pj * ( /• ) ,  making  the  measurement  incorrect,  u 
threshold  was  placed  on  the  output  signal-to-noise 
ratios  of  both  S(k,  i )  and  S'(Ar  —  1,  <).  '1’his,  thresh¬ 
old  required  that  the  -ignal-to-uoise  ratios  on  both 
S(k,  i )  and  ,S’(F— I,  >)  be  greater  than  +  10  dh  if 
p:*(k)  was  to  he  included  in  the  data.  When  this 
rci|iiirement  is  applied,  approximately  100  percent 
of  the  short-range  data,  50  percent  of  the  medium- 
range  data,  and  40  percent  of  the  long-range  data 
meet  or  exceed  the  threshold. 

(U)  Figure  21a  depicts  the  histogram  of  correla¬ 
tion  coefficients  *>.*(/•)  for  the  short-range  data. 
The  mean  correlation  coefficient  is  0.S55,  with  a 
standard  deviation  of  O.IMS.  In  Figs.  21b  and  21c, 
the  corresponding  histograms  for  the  medium-  and 
long-range  data  •ire  shown.  The  mean  correlation 
coefficient  for  lie  medium  data  i'  0.91.1,  with  a 
standard  dev  iation  of  0.112.  The  mean  correlation 
for  the  long-range  data  is  0.91(1  with  a  standard 
deviation  of  0.05-1.  The  slight  increase  in  mean 
correlation  value  for  the  long-range  data  is  not 
considered  to  be  significant, 

(U)  The  dependence  of  spatial  stability  on  the 
relative  direction  of  the  spatial  displacement  is  of 
interest.  The  prevailing  conjecture  is  that  the 
slahilitv  will  he  less  on  a  full  Doppler  track  than 
on  n  /cio  Doppler  track.  Figure  22a  depicts  a 
fy.*ikl  histogram  for  the  long-range  data  on  a 
full  Doppler  track;  Fig.  22b  coricspouds  to  a  zero 
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Figure  22l>.  Histogram  of  pi*(k)  for  zero  Doppler  condi¬ 
tions.  (Figure  uiiclnssitied.) 

Doppler  track.  The  mean  correlation  coefficient 
for  the  zero  Doppler  data  is  0.925,  with  a  standard 
deviation  of  0.051;  for  the  full  Doppler  track  the 
mean  and  standard  deviations  are  0.014  ami  0.05.1, 
respectively.  Thus,  the  January  1974  data  do  not 
indicate  any  dependence  of  spatial  stability  on 
the  direction  of  displacement,  at  least  for  long 
ranges. 

(U)  The  results  of  the  January  1974  spatial 
stability  experiment  are  remarkably  similar  to 
those  obtained  in  the  M5  and  M<>  experiments. 
Figures  2.‘lu  and  214b  compare  the  overall  (short 
medium,  and  long  range)  histogram  of  correlation 
coefficients  obtained  in  the  January  1974  experi¬ 
ment  with  a  typical  one-day  histogram  from  the 
M(K'  experiment.  Although  the  measurement 
interval  was  shorter  in  the  January  experiment 
than  in  the  Mf.C  experiment  (29  see  versus  01  sec ), 
the-  similarities  in  the  histograms  are  amazing, 
considering  the  differences  in  location  and  range 
between  the  two  experiments. 

VI.  FUTURE  APPLICATIONS  AND 
CONCLUSIONS 

A.  CONCLUSIONS 

((.')  The  effectiveness  of  coherent,  matched-filter 
teeliiiicpies  was  demonstrated  in  fixed-site  experi¬ 
ments  in  the  Strait;,  of  Florida.  The  application  of 
coherent  integration  over  time  intervals  of  the 
order  of  one  minute  yielded  satisfactory  operation 
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Fillin'  2:!:i.  CiunpurLon  of  correlation  coefficient  histo. 
Krsmi- :  overall  p3*tfct  histogram  from  January  1974 
experiment.  (Figure  unclassified.) 


til  input  signal-to-noise  ratios  below  0  db.  By 
matched  libeling  of  (lie  received  signal,  tin*  effects 
of  selective  fading  and  intersymhol  interference 
dui'  to  multipath  were  reduced.  Massive  amounts 
of  data  (over  4.000  hours)  on  system  performance 
add  significance  to  the  result-,  obtained. 

((')  To  apply  coherent,  matched-filter  techniques 
to  the  submarine  communication  problem,  an 
understanding  of  the  medium’s  spatial,  as  well  as 
temporal,  stability  is  required.  A  ten-day  experi¬ 
ment  in  the  Atlantic  between  Kleutheru  and 
Bermuda  with  a  towed  source  and  a  fixed  hydro¬ 
phone  ua-  conducted  to  measure  spatial  stability. 
The  results  from  this  experiment  indicated  the 
presence*  of  sufficient  spatial  stability  over  inter¬ 
vals  of  at  least  -SO  sec  tit  six  knots,  at  ranges  from 
0  to  400  nmi. 

((')  The  combination  of  fixed-site  communication 
system  results  and  the  spatial  stability  measure¬ 
ments  establishes  the  feasibility  of  coherent, 
matched-lilter  techniques  for  submarine  communi¬ 
cations.  Advantages  to  be  gained  from  these 
technique'  include  reliable  operation  at  low 
signal-to-noise  ratios  and  under  varying  propaga¬ 
tion  conditions.  Further,  randomi/.cd  transmission 
formats  can  be  employed  to  reduce  the  detectabil¬ 
ity  of  the  communication  signal.  Because  the 
technology  required  to  implement  coherent, 
matched-lilter  techniques  is  readily  available, 
their  application  to  submarine  communications 
should  be  initiated. 
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Figure  231).  Comparison  of  correlation  coefficient  histo¬ 
grams:  one-day  p(k)  histogram  from  MOC  experiment. 
(Figure  unclassified.) 

1.  Subsidiary  Conclusions 

(U)  'Hu*  experimental  program  yielded  several 
secondary  results  that  should  be  considered  in 
subsequent  submarine  communication  system 
design. 

(C)  First,  the  acoustic  medium  exhibits  signifi¬ 
cant  variance  in  input  signal-to-noise  ratio  SNR|„, 
even  under  fixed-site  conditions.  Standard  devia¬ 
tions  of  SNR,,,  of  the  order  of  5  <lb  can  generally 
be  expected.  These  variations  require  that  any 
reliable  communication  system  be  able  to  operate 
successfully  over  a  wide  range  of  input  signal-to- 
noise  ratios.  Coherent,  matched-filter  systems 
satisfy  this  requirement  more  closely  than  the 
incoherent  systems  currently  in  use. 

(C)  Second,  impulsive  noise  must  be  accounted 
for  in  any  system  in  which  low  (PE  <0-001)  bit 
error  probabilities  are  required.  The  presence  of 
infrequent,  high-energy  noise  pulses  can  limit 
the  error  probability  to  a  fixed  level,  independent 
of  the  average  input  signal-to-noise  ratio.  Soft 
limiting  of  the  receiver  input  in  conjunction  with 
error-correcting  codes  cun  be  used  to  overcome 
this  limitation. 

H.  SUGGESTED  APPLICATION 

(C)  The  M5  and  M6  communication  systems 
were  designed  specifically  for  the  evaluation  of 
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Tahi.e  A.  I'.ipiclnl  /u rformuiict  fur  the  M7  .system; 
l\  i  T  In  msec,  U  -  lull  Itz.  ( Table 

(inssijitil  ( 'onjiili ntial.) 


K 

s 

s 

Race , 
hits/sec 

Minimum 
Operating 
SNR,  db 

Range, 

nml 

too 

40 

1.0 

-  4 

<100 

200 

20 

0.5 

-  7 

MOO 

.00 

10 

0.25 

-10 

>100 

I  100 

4 

0.1 

-14 

>100 

coherent ,  matched-filter  techniques.  Consequently, 
their  structures  luck  many  features  required  of  a 
priicl i(*ii !  submarine  communication  system.  For 
<\\um|>li\  tin-  transmission  in  tin-  M5  and  M6 
.systems  is  continuous,  and  synchronization  is 
performed  by  the  operator.  Neither  of  these 
characteristics  is  acceptable  in  practice.  The  pur¬ 
pose  of  this  section  is  to  suggest  one  realization  of 
a  coherent,  inatched-lilter  communications  system 
that  is  suitable  for  operational  use.  Only  the 
structure  of  the  transmission  and  the  expected 
performance  characteristics  are  given  here.  The 
details  of  the  receiver  processing  and  implementa¬ 
tion  will  he  available  in  the  future. 

((’)  Although  good  performance  at  low  signal-to- 
noisc  ratios  and  under  varying  multipath  condi¬ 
tions  is  necessary  for  any  submarine  communica¬ 
tion  s\stem,  other  operationally  important  re¬ 
quirements  must  also  be  considered.  Two  such 
requirements  are  that  the  transmission  be  of 
limited  duration  and  that  it  be  detection  resistant. 
fI  he  first  requirement  precludes  blanking  of 
passive  sonar  systems  on  the  transmitting  plat¬ 
form,  while  the  second  reflects  the  need  to  prevent 
detection  or  localization  of  the  transmitting 
platform. 

((  )  The  M7  system,  des-ribed  below,  has  been 
designed  to  provide  an  application  of  coherent, 
matched-filter  techniques  that  satisfy  the  require¬ 
ments  previously  staled.  Specifically,  the  M7 
system  utilizes  a  hurst  type  of  transmission  with  a 
randomized  lormai.  The  structure  of  the  trans¬ 
mission  is  based  entirely  on  an  arbitrary  random 
binary  sequence  that  can  be  changed  on  a  daily 
basis  il  desired  The  only  constraint  on  this 
sequence  is  that  it  must  he  known  to  both  tin* 


transmitter  and  receiver  and  have  the  same  statis¬ 
tical  properties  as  other  random  sequences. 

(C)  in  the  M7  system,  the  probe  and  information 
components  are  phase  multiplexed  instead  of  time 
multiplexed  as  in  the  Mfl  and  Mb  systems.  Lei 
</(/)  he  a  low -pass  digit  of  duration  7’ and  let  f„  he 
the  carrier  frequency,  then  the  M7  transmission 
is  given  by 

Y}  qr(i)u(t—iT)  cos  2 itfj 

1^0 

-f  f/, (<)«(/ -  ?T;  sin  2r  fut  (So) 

The  first  term  in  the  sum  is  the  probe  component 
and  the  second  is  the  information  component.  The 
coefficients  (/,.:<).  i/i('U|±i)  arc  derived  from  a 
random  binary  sequence  </(!)  and  the.  information 
component  us  indicated  below, 

(C)  Let  </(/)« J  f  I  j  lx-  a  random  binary  sequence 
of  length  2Kb  'Die  coefficients  for  tin*  probe 
component «/,.( >)  are  simply  the  /(<)'sof  even  index; 

'jr(')  =-•</(  2/) ;  A’ -I  (80) 

Assume  that  A',,  information  hits  e\i)  are  to  he 
transmitted.  Further,  assume  that  Ab  is  a  factor 
of  A i  so  tlint  K,  =  K„  Kn  for  some  integer  A',. 
Then  the  coefficients  for  the  information  com¬ 
ponent  are  given  bv 

,m;)=,(|;7yb)h(2;-h);  ;=o  . . .  K,~  i  us:; 

where  ls\  it. dilates  the  integer  part  of  i.  'l  ints, 
7i (  l  is  formed  by  the  inoduhition  of  </(/)  for  odd 
indie:  s  by  the  information  values  in  groups  of  A", 
digits. 

(C)  Since  Lot li  coefficients  (/,.(/)  and  ijt{i)  are 
derived  from  a  random  sequence,  </(<!,  the  trans¬ 
mission  m(t)  can  be  expected  to  be  free  from 
periodicities  or  strong  spectral  components,  >f 
the  input  signal-to-noise  rntio  of  an  intercept 
receiver  is  below  0  db,  the  detectability  of  the  M7 
transmission  will  he  verv  low.  if  tie*  input  signttl- 
to-nois-e  ratio  of  the  intercept  receiver  is  above  U 
db,  eonv ‘iitional  power  measurement  techniques 
are  applicable.  ( 'onscquontly,  the  M7  leehnicpie 
(or  any  other)  can  he  considered  detection  resistant 
only  if  the  intercept  receiver  is  denied  an  adequate 
signal-to-noi'C  ratio,  that  is,  he  beyond  a  certain 
rauge. 
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(C)  To  provide  nn  Illustration  of  tin*  performance 
to  l«'  i*.\ peel i‘d  of  tin-  M7  system,  a  specific 
practical  example  is  given  Assume  that  the  system 
bandwidth  is  10(1  hz  and  that  the  transmission 
duration  is  limited  to  40  see.  '.cl  </!/)  be  a  10- 
msec  reel  angular  pulse1,  so  that  the  spectrum  of 
//(()  (ills  tlu*  system  bandwidth.  Then,  K <  espials 
4,000.  Table  5  gives  the  system  performance  for 
several  choices  of  A’,.  The  column  marked  “mini¬ 
mum  operating  SNR”  gives  the  input  signal-to- 
noise  ratio  expected  to  yield  a  bit  probability 
of  error  of  0.001.  This  signal-to-noise  ratio  is 
obtained  by  taking  the  theoretical  signal-to- 
noise  ratio  required  to  obtain  such  performance 
and  then  adding  a  +0  db  differential  to  account 
for  nonidealities.  Note  the  low  signal-to-noise 
ratios  for  which  the  0.001  error  probability  is 
expected. 
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